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Abstract 

The AcerMC Monte Carlo Event Generator is dedicated for the generation of Standard Model background processes 
at pp LHC collisions. The program itself provides a library of the massive matrix elements and phase space modules for 
generation of selected processes: gg,qq tibb; qqW{-^ £i/)bb; qqW{~^ lv)tt; gg,qq — * ^/7*(— > U)bb; gg,qq —> 
Z/'y*{^ U,viy,bb)tt; complete electroweak gg,qq {Z/W/y* —^)bbtt; gg,qq ttti; gg,qq {tt —*)ffbffb; 
gg, qq (WWbb ~*)ffffbb. The hard process event, generated with one of these modules, can be completed by the initial 
and final state radiation, hadronisation and decays, simulated with either PYTHIA or HERWIG Monte Carlo event generator 
and (optionally) with TAUOLA and PHOTOS. Interfaces to all these packages are provided in the distribution version. The 
matrix element codes have been derived with the help of the MADGRAPH package. The phase-space generation is based on 
the multi-channel self-optimising approach using the modified Kajantie-Byckling formalism for phase space construction and 
further smoothing of the phase space was obtained by using a modified ac-VEGAS algorithm. 
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1 PROGRAM SUMMARY 

Title of the program: AcerMC version 2.0 

Operating system: Linux 

Programming language: FORTRAN 77 witli popular extensions. 
External libraries: CERNLIB. 

Size of the compressed distribution directory: about 10.7 MB. The distribution includes modified versions of PYTHIA 6.2, 
HERWIG 6.5, HELAS libraries, TAUOLA and PHOTOS packages. 

Key words: Standard Model backgrounds at LHC, massive matrix elements, Monte Carlo generator, heavy flavor production, 
multi-channel phase-space generation. 



Nature of physical problem: Despite a large repertoire of processes implemented for generation in event generators like 
PYTHIA [1] or HERWIG [2] a number of background processes, crucial for studying expected physics potential of the LHC 
experiments is missing. For some of these processes the matrix element expressions are rather lengthly and/or to achieve a 
reasonable generation efficiency it is necessary to tailor the phase-space selection procedure to the dynamics of the process. 
That is why it is not practical to imagine that any of the above general purpose generators will contain every, or even only 
observable, processes which will occur at LHC collisions. A more practical solution can be found in a library of dedicated 
matrix-element-based generators, with the standardised interfaces like that proposed in [3], to the more universal one which 
is used to complete the event generation. 

Method of solution: The AcerMC Event Generator provides itself library of the matrix-element-based generators for several 
processes. The initial- and final- state showers, beam remnants and underlying events, fragmentation and remaining decays 
are supposed to be performed by the other universal generator to which this one is interfaced. We will call it supervising 
generator. The interfaces to PYTHIA 6.2 and HERWIG 6.5, as such generators, are provided. Provided is also interface to 
TAUOLA [4] and PHOTOS [5] packages for r-lepton decays (including spin correlations treatement) and QED radiations 
in decays of particles. At present, the following matrix-element-based processes have been implemented: gg,qq ttbb, 
qq W{-^ £iy)bb; qq W{-^ tu)tt\ gg,qq Z/-y*{-^ U)bb\ gg,qq Z/f*{-^ ££,uiy,bb)ti- complete EW 
gg,qq ^ {Z/W/^* ^)tibb- gg,qq ^ titi; gg,qq ^ {tt -^)ffbffb- gg,qq ^ {WWbb ^)ffffbb. Both interfaces 
allow the use of the PDFLIB library of parton density functions. Provided is also set of control processes: qq ^ W ^ 
qq Z/j* U\ gg, qq and gg {tt ^)WbWb- 

Restriction on the complexity of the problem: The package is dedicated for the 14 TeV pp collision simulated in the LHC 
environment. In particular, using it for the 2 TeV pp collision although technically possible might not be efficient and would 
require dedicated optimisation of the phase space generation. The consistency between results of the complete generation 
using PYTHIA 6.2 or HERWIG 6.5 interfaces is technically limited by the different approaches taken in both these generators 
for evaluating uqcd and agsD couplings and by the different models for fragmentation/hadronisation. For the consistency 
check, in the AcerMC library contains native coded definitions of the agco and uqed- Using these native definitions leads 
to the same total cross-sections both with PYTHIA 6.2 or HERWIG 6.5 interfaces. 

Typical running time: On an PHI 800 MHz PC it amounts to ~ 0.05 — ^ LI events/sec, depending on the choice of process. 

[1]. T. Sjostrand et al.. High energy physics generation with PYTHIA 6.2, eprint hep-ph/0108264 LU-TP 01-21, August 
2001. 

[2]. G. Marchesini et al., Comp. Phys. Commun. 67 (1992) 465, G. Corcella et al, JHEP 0101 (2001) 010. 
[3]. E. Boos at al.. Generic user process interface for event generators, hep-ph/0 109068 

[4]. S. Jadach, J. H. Kuhn, Z. Was, Comput. Phys. Commun. 64 (1990) 275; M. Jezabek, Z. Was, S. Jadach, J. H. Kuhn, 
Comput. Phys. Commun. 70 (1992) 69; R. Decker, S. Jadach, J. H. Kuhn, Z. Was, Comput. Phys. Commun. 76 (1993) 361. 
[5]. E. Barberio and Z. Was, Comp. Phys. Commun. 79 (1994) 291. 
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2 Changes since AcerMC 1.0 

• AcerMC version 1.1: The changes include transition to HERWIG 6.4, updated scale choices for processes 5-8 (c.f. 
Section ItTtI and the inclusion of control processes 91-94 for consistent process evaluation (c.f. Section|49}. Also, a 
possibility of an event dump according to the Les Houches standard was added (see Section lT!9l . 

• AcerMC version 1.2: A bug fix in HERWIG 6.4, affecting the shower evolution, was made. It carries no immediate 
impact on the AcerMC processes but was discovered by the AcerMC authors and added for the convenience of the users. 
This bug fix will be included in future versions of HERWIG. Also, the PYTHIA version was upgraded to PYTHIA 
6.208 and the the implementation of storing/reading back of hard process events according to the Les Hauches standard 
was simplified with respect to AcerMC 1.1 (see Section FTgl . Also, for the convenience of the users the Pythla code 
was modified so that the top decay products from AcerMC processes are now stored in the history part of the event 
record (status code 21) and have the correct pointers to the top quark they originate from. This feature will be added to 
the future versions of PYTHIA. 

• AcerMC version 1.3: A transition to the HERWIG 6.5 was made, which now supports the Les Houches standard 
for handling the external processes. Consequently, the AcerMC interface to HERWIG was completely rewritten (c.f. 
Section l7!6l . As a direct consequence, the same event record can freely be swapped (i.e. read back) to either PYTHIA or 
HERWIG for fragmentation and hadronisation treatment. In addition, the PYTHIA version was upgraded to PYTHIA 
6.214. Also, the build procedure of the libraries and executables was greatly simplified (see Section lTlsl . 

• AcerMC version 1.4: The interfaces to external TAUOLA and PHOTOS libraries were added. The necessary modifi- 
cations in the interface routines and the native PHOTOS code were made to enable the user to process the events with 
TAUOLA and/or PHOTOS using PYTHIA or HERWIG as the supervising generators. 

• AcerMC version 2.0: New algorithm for phase space generation implemented and optimised. 

New processes were added: qq — > [Z/W/'y* —*)ttbb\ gg,qq — > tttt; gg,qq — > {tt ~*).f fbf fb\ gg,qq — > ti and 
gg {tt ^)WbWb. New control channel added: gg [tt ^)WbWb. 
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3 Introduction 

Despite a large repertoire of processes implemented for generation in the universal generators like PYTHIA (V\ or HERWIG 
|2) a number of Standard Model background processes for studying expected physics potential of the LHC experiments is 
still missing. For some of these processes the matrix element expressions are rather lengthy and/or to achieve a reasonable 
generation efficiency it is necessary to tailor the phase-space selection procedure to the dynamics of the process. Due to this 
fact it cannot be expected that any of the universal purpose generators will contain every proces that is expected to participate 
at LHC pp collisions with an observable rate. A more practical solution could come in form of dedicated matrix-element-based 
generators with standardised interfaces, like the one proposed in to the more general ones which are used to complete 
event generation. 

The AcerMC Monte Carlo Event Generator follows up on this idea. It is dedicated for the simulation of the specific 
Standard Model background processes at LHC collisions: the gg,qq —> tibb, qq W{~^ iv)bb; qq W{~^ iv)ti; 
gg,qq Z/^*{-^ ££)bb; gg,qq ^ Z/j*{-^ U,vv,bb)ti; complete EW gg,qq —> {Z/W/^* ^)ttbb- 4 top-quark 
production gg,qq tttt; extended treatement of the 2 top-quark production gg,qq {ti -^)ffbffb and gg,qq — * 
{WWbb —*)ffffbb. They are characterised by the presence of the heavy flavour jets and multiple isolated leptons in 
the final state. For the Higgs boson searches, the tiH, ZH, WH with H bb, the gg -* H with H ZZ* U, 
the bbh/H/A with h/H/A — > tt,^^ are the most obvious examples of signals where the implemented processes would 
contribute to the dominant irreducible backgrounds. The same background processes should also be considered for e.g. 
estimating the observability of SUSY events with a signature of multi-b-jet and multi-lepton production. 

The program itself provides library of the massive matrix elements and phase space modules for the generation of the 
implemented processes. The hard process event, generated with these modules, can be completed by the initial and final state 
radiation, hadronisation and decays, simulated with either PYTHIA 6.2 or HERWIG 6.5 Monte Carlo Event Generators. 
These will subsequently be called the Supervising Generators. Interfaces of AcerMC to both, PYTHIA 6.2 or HERWIG 6.5 
generators, are provided in the distribution version. Provided is also interface to TAUOLA 1 5 1 and PHOTOS 1 4 1 packages, for 
the more correct treatement of the r-lepton decays and photon radiation, than what available in the Super-vising Generators. 
The AcerMC also uses several other external libraries: CERNLIB, HELAS |6 1, VEGAS |7|. The matrix element codes have 
been derived with the help of MADGRAPH |8| package. The achieved typical efficiency for the generation of unweighted 
events is of 20% - 30%, rather high given a complicated topology of the implemented processes. 

This paper superseeds the first version of the manual, published in |9|. The outline of this paper is as follows. In Section 
3, we describe physics motivation for implementing each of the above processes and we collect some numerical results (plots, 
tables) which can be used as benchmarks. In Section 4 we describe the overall Monte Carlo algorithm. Section 5 gives 
details on the structure of the program. Section 6 collects information on how to use this program and existing interfaces to 
PYTHIA 6.2 and HERWIG 6.5, TAUOLA and PHOTOS. Summary, Section 7, closes the paper. Appendix A documents 
sets of Feynman diagrams used for calculation of the matrix element for each subprocess. Appendices B and C give examples 
of the input/output of the program. 
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4 Physics content 

The physics programme of the general purpose LHC experiments, ATLAS [14] and CMS (15], focuses on the searches for 
the New Physics with the distinctive signatures indicating production of the Higgs boson, SUSY particles, exotic particles, 
etc. The expected environment will in most cases be very difficult, with the signal to background ratio being quite low, on the 
level of a few percent after final selection in the signal window. 

Efficient and reliable Monte Carlo generators, which allow one to understand and predict background contributions, are 
becoming the key point to the discovery. As the cross-section for signal events is rather low, even rare Standard Model 
processes might become the overwhelming background in such searches. In several cases, generation of such a process 
is not implemented in the general purpose Monte Carlo generators, when the complicated phase space behaviour requires 
dedicated (and often rather complex) pre-sampling, whilst the general purpose Monte Carlo generators due to a large number 
of implemented processes tend to use simpler (albeit more generic) phase space sampling algorithms. In addition, the matrix 
element for these processes is often lengthy and thus requiring complicated calculations. Only recently, with the appearance 
of modem techniques for automatic computations, their availability on demand became feasible for the tree-type processes. 
With the computation power becoming more and more easily available, even very complicated formulas can now be calculated 
within a reasonable time frame. 

The physics processes implemented in AcerMC library represent such a set of cases. They are all being key background 
processes for the discovery in the chaimels characterised by the presence of the heavy flavour jets and/or multiple isolated 
leptons. For the Higgs boson searches, the ttH, ZH, WH with H bb, the gg ^ H with H ZZ* ~+ 4^, the bbh/H/A 
with h/H/A TT, fifi are the most obvious examples of such channels. 

It is not always the case that the matrix element calculations in the lowest order for a given topology represent the total 
expected background of a given type. This particularly concerns the heavy flavour content of the event. The heavy flavour 
in a given event might occur in the hard process of a much simpler topology, as the effect of including higher order QCD 
corrections (eg. in the shower mechanism). This is the case for the b-quarks present in the inclusive Z-boson or W-boson 
production, which has a total cross-section orders of magnitude higher than the discussed matrix-element-based Wbb or Zbb 
production. Nevertheless, the matrix-element-based calculation is a very good reference point to compare with parton shower 
approaches in different fragmentation/hadronisation models. It also helps to study matching procedures between calculations 
in a fixed agco order and parton shower approaches. For very exclusive hard topologies matrix-element-based calculations 
represent a much more conservative approximation than the parton shower ones 1 16 1. 

Let us shortly discuss the motivation for these few Standard Model background processes which are implemented in the 
AcerMC 2.0 library. 

The ttbb production at LHC is a dominant irreducible background for the Standard Model (SM) and Minimal Super- 
symmetric Standard Model (MSSM) Higgs boson search in the associated production, tiH, followed by the decay H — > bb. 
The potential for the observability of this channel has been carefully studied and documented in 1 14| and 1 17|. Proposed 
analysis requires identifying four b-jets, reconstruction of both top-quarks in the hadronic and leptonic mode and visibility of 
the peak in the invariant mass distribution of the remaining b-jets. The irreducible ttbb background contributes about 60-70% 
of the total background from the ti events {ttbb, tibj, tijj). 

The Wbb production at LHC is recognised as a substantial irreducible background for the Standard Model (SM) and 
Minimal Supersymmetric Standard Model (MSSM) Higgs boson search in the associated production, WH, followed by the 
decay H — > bb. The AcerMC library discussed here includes even more efficient implementation of the algorithm presented 
in(l^. 

The Wtt production at LHC has to our knowledge not been implemented in the publicly available code so far' . It 
is of interest because it contributes an overwhelming background 1 13| for the measurement of the Standard Model Higgs 
self-couplings at LHC in the most promising channel pp HH WWWW . 

The Z — > £l)bb production at LHC has since several years been recognised as one of the most substantial ir- 
reducible (or reducible) backgrounds for the several Standard Model (SM) and Minimal Supersymmetric Standard Model 
(MSSM) Higgs boson decay modes as well as for observability of the SUSY particles. There is a rather wide spectrum of 
regions of interest for this background. In all cases the leptonic Z/'y* decay is asked for, but events with di-lepton invariant 
mass around the mass of the Z-boson mass or with the masses above or below the resonance peak could be of interest. The 
presented process enters an analysis either by the accompanying b-quarks being tagged as b-jets, or by the presence of leptons 
from the b-quark semi-leptonic decays in these events, in both cases thus contributing to the respective backgrounds. 

Good understanding of this background, and having a credible Monte Carlo generator available, which allows studying 
of expected acceptances for different final states topologies, is crucial for several analyses at LHC. The AcerMC library 

'We thank M. L. Mangano for bringing this process to our attention and for providing benchmark numbers for verifying the total cross-section. 
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discussed here includes more efficient implementation of the algorithm presented in fllV 

The Z 1^* ( — > it, vv, bb)tt production at LHC is an irreducible background to the Higgs search in the invisible decay 
mode (case of Z ^ vu) in the production with association to the top-quark pair |28|. With the Z/'y*[~* bb) it is also an 
irreducible resonant background to the Higgs search in the ttH production channel but with the Higgs boson decaying to the 
b-quark pair 1 17 1. 

The complete EW production of the gg, qq — > {Z/W/'^* — >)b6tt final state is also provided. It can be considered 
as a benchmark for the previous process, where only the diagrams with resonant gg,qq —> {Z/^* — »)b6tt are included. It 
thus allows the verification of the question, whether the EW resonant contribution is sufficient in case of studying the ttbb 
background away from the Z-boson peak, like for the ttH with Higgs-boson mass of 120 GeV. 

The gg, qq — »■ tttt production , interesting process per se, is a background to the possible Higgs self-coupling mea- 
surement in the gg HH WWWW decay, 1 13 1. 

Thegg,qq — > {WWbb — >) f f f fbb and gg, qq — > {tt — processes give posiblity to study spin 
correlations in the top-quark pair production and decays as well as the effect from the off-shell production. Those are important 
for the selection optimisation eg. in the gg H ^ WW channel, see the discussion in f351 

A set of control channels, i.e. the qq —>■ Z/~f* II, gg, qq ti , qq ^ W ^ iv and gg — > (ti -^^WhWb 
processes, have been added to AcerMC in order to provide a means of consistency and cross-check studies. 

This completes the list of the native AcerMC processes implemented so far. Having all these different production pro- 
cesses implemented in the consistent framework, which can also be directly used for generating standard subprocesses imple- 
mented in either PYTHIA or HERWIG Monte Carlo, represents a very convenient environment for several phenomenological 
studies dedicated to the LHC physics. 

For the cases, where radiative photon emission from final state leptons is important the package PHOTOS O can be used 
in the chain of event generation. In similar way also package TAUOLA |5 1 can be interfaced directly to the generation chain 
and used for events generation in cases where more detailed treatment of the tau-lepton decay and including spin correlations 
effects is relevant. 

In the following subsections we discuss in more detail implementation of each subprocess. We also give benchmark Tables 
with the total cross-sections obtained with AcerMC processes but different implementations and setting of qqcd(Qqcd): 
the native AcerMC, PYTHIA and HERWIG ones. For a more detailed discussion on this topic the reader is referred to Section 
15.6! If the native AcerMC definition is used, the same cross-section is obtained either with PYTHIA or HERWIG generation 
chains. 

4.1 The gg,qq ^ tibb processes 

In the implementation discussed here, the matrix element was derived using the MADGRAPH package (s). These matrix 
elements are not covering the decay of the top-quarks, the latter are considered as massive final states of the process. The 
top-quark decays is than performed by the supervising generator. Rather important spin effects (spin correlations) in the top 
decays are therefore not yet included. The similar solution, like for tau decay in the Z-boson production process discussed in 
1291 . is planned to be implemented here in the near future. 

As a benchmark, the processes gg, qq ttbb have been simulated for pp collisions with 14 TeV centre-of-mass energy 
and CTEQ5L 1 30 1 parton density functions, using event generation with massive 2-^4 matrix element implemented as an 
external process to PYTHIA 6.2 (see Section 4 and 5). The decays of the top-quarks have been left under control of PYTHIA 

6.2 generator. The qq ttbb subprocess contributes less than 10% of the total cross-section. 

The total cross-section is very sensitive to the choice of the QCD energy scale used for calculation of that process, thus 
indicating potentially large contributions from higher order corrections. The same definition for the factorisation and renor- 
malisation scale is used. The example values of the total cross-section for implemented choices of the QCD energy scale are 
given in TableEf. 

As a cross-check, the processes gg, qq — > ttbb have been coded independently using the COMPHEP package |31 1. The 
same set of diagrams was selected and only the integrating part of the package was used to calculate total cross-section. The 
choices for the QCD energy scale were kept consistent. A very good agreement between the cross-sections obtained with two 
independent calculation streams prepared for this study has been achieved II2I . 

One can observe a very strong scale dependence of the cross-section for the gg, qq — > ttbb process (c.f. Table0. Factor 
four (! !) can be expected on the predicted cross-section when changing from the scale Qqcd = s to the scale Qqqd ~ < 

^Numbers obtained with HERWIG generator (third collumn) are slightly different than what published in jO). This is related to the internal changes in 
HERWIG between version 6.3 and 6.5. 
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Table 1; AcerMC cross-sections for the gg, qq — > tibb production at different choices of the 
QCD energy scale and aqcD implementations. The 14 TeV centre-of-mass energy 
and CTEQ5L parton density functions were used for the simulation with interfaces to 
PYTH\A 6.2 and HER\N\G 6.5. The ruH = 120 GeV and nit = 175 GeV were used 
for calculating the Qqqu in the last row of this table. The default settings of uqcd 
as implemented in AcerMC, PYTHIA 6.2 and HERWIG 6.5 were used. 



Factorisation scale 


OIQCD (IL) 
native AcerMC 


OQCD (IL) 

as in PYTHIA 6.2 


O-QCD (2L) 

as in HERWIG 6.5 




gg — * tthh 


Qqcd = ^ 


4.2 [pb] 


3.9 [pb] 


3.0 [pb] 


Qqcd = E (Pt' + ™?)/4 


10.3 [pb] 


10.2 [pb] 


7.2 [pb] 


QlcD = E ipV)/i 


17.0 [pb] 


16.9 [pb] 


11.5 [pb] 


Qqcd = (™t + ™h/2)^ 


8.2 [pb] 


8.1 [pb] 


5.8 [pb] 




qq — > tthh 


Qqcd = ^ 


0.30 [pb] 


0.29 [pb] 


0.22 [pb] 


Qqcd = T. {PT" + "^f}/'i 


0.61 [pb] 


0.60 [pb] 


0.43 [pb] 


QlcD = Y. (Pt')/4 


0.91 [pb] 


0.90 [pb] 


0.62 [pb] 


Qqcd = (™t + ™h/2)^ 


0.52 [pb] 


0.51 [pb] 


0.37 [pb] 



>. This very strong dependence on the energy scale is also observed in the case of the ttH production, for recent 
discussion see J32\. There, the recommended central factorisation and renormalisation energy scale is = {mt + mn/S). 
Having in mind that the primary interest of evaluating this background is the Higgs search in the ttH production, i.e. with 
the b-quark system being produced with the invariant mass of the expected Higgs boson, we have also introduced this central 
energy scale, with tuh ~ 120 GeV as one of the possible choices. 

Fig. Q] shows the distributions of the Qqcd = \JQqcd (distributions have been normalised to one) for the ttbb 

events with the invariant mass of the b-jets system, calculated using the default PYTHIA (LO) aqcD implementation^. 
rribb-jets = 120 ± 30 GeV. The distribution is well collimated around the average value when Qqcd is defined as 
< my > or < p|n > while it is much broader when Qqcd is defined as s. The kinematic distributions are very similar in 
shape for separate gg ttbb and qg -+ tibb contributions. The total cross-section for accepted events as a function of the 
averaged Qqcd (for these events) is shown in the bottom plot. It can be noted that the cross-section decreases rather fast with 
the increasing value of the average < Qqcd >■ Also shown is the ct^iQQCD) dependence scaled to match the cross-section 
at Qqcd = (rnul'i + vat) with mu = 120 GeV, it being the only calculated cross-section point with a fixed scale. The 
other cross-sections are shown to follow the expected at (Qqcd) dependence rather well, while the deviations are induced 
by the parton density function dependence on the Qqcd scale, most notably at Qqcd = ^ value. The deviations induced by 
the parton density functions dependence on the Qqcd scale are different for the gg and qq contributions, as can be concluded 
from results given in TablelB 

The series of plots illustrating the most relevant differential distributions for the top-quarks and b-quarks can be found in 

4.2 The qq ^ lv)g*{-^ hh) process 

The matrix element for the implemented process was coded by using the MADGRAPH package fsl. This process is repre- 
sented by only two Feynman diagrams, with quark exchange in the t-channel, leading to the production of the VF-boson and 
virtual gluon splitting into bb pair. Only the u, d, s, c quarks were considered in this implementation, the possibility of the 
b-quark in the initial state was omitted as expected to be negligible numerically (e.g.. \Vbc/Vud\^ ^ 0.002) but leading to 
several additional diagrams which would have to be included. The massive matrix element takes into account spin correlations 
in the W-boson decay and angular correlations between leptons and quarks. Due to the massive treatment of the final state 

^This would makes distributions directly relevant for the ttH analysis. For details on the jet reconstruction see |33 . 
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Fig. 1: Top: the Qqcd distributions for ttbb events with the invariant mass of the b-jets system 
fTibb-jets = 120 ± 30 GeV. Bottom: the total cross-section of accepted events as a 
function of the averaged Qqcq (for these events). 



fermions the amplitude has no singularities; the total cross-section is well defined. The effect from the -boson natural width 
and the -boson propagator are also properly included. 



Table 2: AcerMC production cross-sections for the qq — > Whh with W ^ ev decay (single 
flavour). The 14 TeV centre-of-mass energy and CTEQ5L parton density functions 
were used with different definitions of oiqed, olqcd ( i^s in default PYTHIA 6.2 and 
HERWIG 6.5J and several choices of the factorisation scale, aqED and aqcD imple- 
mentations. 



Factorisation scale 


OiQED, O-QCD (IL) 

native AcerMC 


O-QED, OLQCD (IL) 

as in PYTHIA 6.2 


aQED, OLQCn (2L) 

as in HERWIG 6.5 




36.5 [pb] 


36.4 [pb] 


30.6 [pb] 




44.1 [pb] 


44.0 [pb] 


36.0 [pb] 




36.0 [pb] 


36.0 [pb] 


29.8 [pb] 


= {s^y + s:s)/2 + pT^y 


37.2 [pb] 


37.1 [pb] 


30.4 [pb] 



As a benchmark, the process qq ^ W(-^ (■v)g'' {-^ bb) has been simulated for pp collision with 14 TeV centre-of-mass 
energy. The total cross-section, including branching ratio for W ^ iv (single flavour) is 36.5 pb (CTEQ5L parton density 
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functions, = Af^„ PYTHIA 6.2 interface)". 

The dependence on the choice of the factorisation scale is rather modest (c.f. Table|2j and does not exceed 20% for the 
choices implemented in AcerMC library. The variation of the cross-section due to different oqed and agco implementa- 
tions and default settings is again evident; as one can expect the two-loop agco implementation given in HERWIG gives a 
~20 % lower cross-section when compared to the cases when native AcerMC and PYTHIA one-loop oiqcd were used^. 

The differential distributions of the qq Wbb events turn out to be interesting when compared to the corresponding ones 
of the qq Zbb and gg — > Zbb events (generated with pure Z-boson exchange). Such comparison is well documented in 

m 

4.3 The qq ^W{^ lv)g*{-^ tt) process 

The 2^4 matrix elements, coded by the MADGRAPH package J8j|, are not covering the decay of the top-quarks; the 
latter are considered as massive final states of the process. The top decay is than performed by the supervising generator. As 
in the case of gg, qq — > tibb process spin effects in the top decays are therefore not yet included. This process, although 
rare, contributes an overwhelming irreducible background to possible measurement of the Higgs-boson self-coupling in the 
HH ^ WWWW decay mode [HI. 

Table 3: AcerMC production cross-sections for the qq — > Wtt with primary W ^ ev decay 
(single flavour). The 14 TeV centre-of-mass energy, CTEQ5L parton density functions 
with different factorisation scales and different deflnitions of the uqed and aqcD 
were used in the matrix element calculations. 



Factorisation scale 


OiQED, OLQCD (IL) 

native AcerMC 


O-QED, OiQCD (IL) 

as in PYTHIA 6.2 


OiQED, OIQCD (2L) 

as in HERWIG 6.5 


QqcD — ^'^tv 


69.3 [fb] 


69. 1 [fb] 


57.4 [fb] 


Qqcd — ^tt 


40.9 [fb] 


39.9 [fb] 


34.7 [fb] 


Qqcd — ^'^w P'^w 


59.7 [fb] 


59.5 [fb] 


49.6 [fb] 


Qqcd = i^w + <t)/'^ + pT'w 


43.7[ft] 


42.8 [fb] 


36.9 [fb] 



Table|3|shows the expected AcerMC cross-sections for different choices of the energy scale and coupling (agED, agco) 
definitions. One should notice the effect of almost a factor two from different choices of the energy scale. 

4.4 The gg, qq — s- Z/7*(^ ££)bb processes 

The matrix elements, derived using the MADGRAPH package 1 8 1, properly take into account spin correlations in the Z-boson 
decay and angular correlations between leptons and quarks. Thank to keeping non-zero b-quark masses the amplitude has no 
singularities; the total cross-section is well defined. 

The full Z/^* exchange proves to be important: For events well below the Z-boson resonance the contribution from 7* 
becomes dominant; the 7* contribution is also sizeable in the high mass tail and increases proportionally with the effective 
mass of the di-lepton system. 

As a benchmark result, the process has been simulated for pp collisions at 14 TeV centre-of-mass energy. The total cross- 
sections, including the branching ratio for Z/j* U (single flavour) are given in Table|4|for different definitions of agsD, 
otgcD couplings. 

Several differential benchmark distributions for leptons and b-quarks originating from the hard process has been collected 
and discussed in 1101 . 



''This can be compared with the matrix element implementation to HERWIG 5.6, u.sed in I34I . I36I . where originally this cross-section was estimated to 
19.8 pb (CTEQ2L parton density functions) but, when implementing CTEQ5L parton density functions and setting kinematic parameters to be in approximate 
accordance with PYTHIA defaults, rises to 36.0 pb, which is consistent with the AcerMC implementation by taldng into account the remaining differences 
in the two calculations (e.g. the former implementation uses an on-shell W boson in the IVIE calculation). 

^ While performing further comparisons of native AcerMC and PYTHIA processes we discovered a misinterpretation of our CKM matrix implementa- 
tion. This connection efectively changes the cross-section for qq — > Wbb and qq — > Wtt processes by '^10% compared to the draft versions of this paper, 
which is nevertheless still well within the physics precision of the program. The affected tables in this paper are already updated. 
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Table 4; AcerMC production cross-sections for the gg, qq — > Z/"f*bb with Z/"f* ee decay 
(single flavour). The 14 TeV centre-of-mass energy, CTEQ5L parton density functions 
and different definitions for the uqed, ctQCD (os in default PYTHIA 6.2 and HERWIG 
6.5j were used in the matrix element calculations. The threshold mu > 10 GeV was 
used in the event generation. 



Factorisation scale 


aQED, aQCD{lL) 
native AcerMC 


aQED, QQcr>(li) 
as in PYTHIA 6.2 


aQED, aQCDC^L) 

as in HERWIG 6.5 




gg Z/j'bb 


= 


49.5 [pb] 


45.8 [pb] 


38.0 [pb] 


= 


53.8 [pb] 


53.9 [pb] 


44.0 [pb] 




54.7 [pb] 


54.6 [pb] 


44.4 [pb] 


Q- = (pT^+s*,)/2 


49.7 [pb] 


49.7 [pb] 


40.5 [pb] 




qq Z/j'bb 




6.7 [pb] 


6.7 [pb] 


5.6 [pb] 


= <S 


8.0 [pb] 


8.0 [pb] 


6.4 [pb] 




7.0 [pb] 


7.0 [pb] 


5.7 [pb] 


Q" = {pT'z + <s)/2 


6.9 [pb] 


6.9 [pb] 


5.7 [pb] 



4.5 The gg, qq — s- Z/7*(^ vv, bb)tt processes 

This process, in spite of having a very small cross-section at LHC energies, contributes as irreducible background to the ttH 
production at low masses. In case the Higgs boson is searched within the H ^ bb mode, this contribution becomes less 
and less important with the Higgs boson mass moving away from the Z-boson mass. In case of the Higgs-boson search in 
the invisible decaying mode, the Z vv might be more relevant also for the higher masses, as the mass peak cannot be 
reconstructed for signal events. The Zj^* U decay is of less interest, as the expected observability at LHC is very low 
(Tablell}. 

Table 5: AcerMC production cross-sections for the gg,qq Ztt with Z —> Vgi^e decay (3 
flavours). The 14 TeV centre-of-mass energy, CTEQ5L parton density functions and 
different definitions for the olqed, ctgcD (os in native AcerMC, default PYTHIA 6.2 
and HERWIG 6.5j were used. 



Factorisation scale 


aQED, aQCDCi-lj) 
native AcerMC 


aQED, aQCD(ii^) 

as in PYTHIA 6.2 


aQED, aQCD['2L) 

as in HERWIG 6.5 




gg I^el^e)it 




126.0 [fb] 


125.8 [fb] 


104. [fb] 




61.9 [fb] 


60.3 [fb] 


52.4 [fb] 




126.2 [fb] 


126.1 [fb] 


105. [fb] 


= pT'i + s,*f)/2 


67.6 [fb] 


66.3 [fb] 


57.0 [fb] 




qq — > Z[ — > iycl^tz)ii 


Q'' = mi. 


64.7 [fb] 


64.6 [fb] 


53.7 [fb] 


= 


39.2 [fb] 


38.2 [fb] 


33.3 [fb] 


Q" = s*z 


64.8 [fb] 


64.7 [fb] 


53.5 [fb] 




41.7 [fb] 


41.0 [fb] 


35.2 [fb] 
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4.6 The electroweak gg, qq (Z/W/^* -^)bbtt process 

One should be well aware, that the gg,qq —> Z/j*tt with Z/"/* — > bb does not represent a complete electroweak production 
of the tibb final state. Consequently, a separate implementation for generation of the complete set of such diagrams (including 
as well W-boson exchange) was addressed. In fact this final state leads to complicated pattern of the 72 Feynman diagrams 
(in case of the gg initial state). 

The contribution from all non-resonant channels is a dominant one for the inclusive cross-section, see Table|S| An almost 
factor 10 higher cross-section is calculated with the full electroweak gg — ^ {Z/W/'y* —*)bbtt with respect to calculated 
with the gg {Z/'y* — > bb)ti process only. One should also note that the electroweak gg — > {Z/W/'y* ~>)bbtt inclusive 
cross-section is on the level of 10% of the QCD gg —> bbtt cross-section, see TableQ] for the same choice of the energy scale. 
But in the mass range around 120 GeV it is on the level of 50% of the QCD contribution, jg). 



Table 6: AcerMC production cross-sections for the electroweak gg, qq (Z/W/^* ~^)bbtt. 

The 14 TeV centre-of-mass energy and CTEQ5L parton density functions were used 
along with different definitions for the aqED, ctQCD (os in native AcerMC, default 
PYTHIA 6.2 fl«t/ HERWIG 6. 5j. The mu = 120 GeV was used for calculation of the 
energy scale. 



Factorisation scale 


OiQED, O-QCD (l-'j) 

native AcerMC 


CQED, aQCD(ii^) 

as inPYTHIA6.2 


oiQED, aQCD{2L) 

as in HERWIG 6.5 




gg ^ {Z/Wh' -^)bbtt 


Qqcd = ^ 


0.58 [pb] 


0.56 [pb] 


0.50 [pb] 


QlcD = Y. {PT" + "^f}/i 


1.10 [pb] 


1.05 [pb] 


0.84 [pb] 


Ql.cn = E {pV)I^ 


1.50 [pb] 


1.50 [pb] 


1.16 [pb] 


Qqcd = (mt + mff/a)^ 


0.90 [pb] 


0.89 [pb] 


0.71 [pb] 




qq (Z/Wh' ~,)bbtt 


Qqcd = ^ 


0.029 [pb] 


0.029 [pb] 


0.025 [pb] 


= I:(Pt'+'"?)/4 


0.043 [pb] 


0.042 [pb] 


0.036 [pb] 




0.049 [pb] 


0.048 [pb] 


0.040 [pb] 


= (mt +m„/2)^ 


0.041 [pb] 


0.041 [pb] 


0.035 [pb] 



4.7 The gg, qq (WWbb -^)ffffbb; gg, qq {tt -^)ffffbb processes 

The implemented 2 — > 6 matrix elements for the resonant gg,qq {tt -^)ffffbb and complete gg,qq —> {WWbb — > 
)ffffbb processes give possibility to study background from top-quark production in more details, than with resonant on- 
shell production only ( as implemented in PYTHIA and HERWIG). In particular, for the Higgs boson search in the H — > 
WW — > Ivif decay channel, the analysis foresees strong suppression against tt background using topological feaures of 
events (jet veto, lepton angluar correlations), but does not foresees implicit top-quarks reconstruction. To reliably predict such 
backgrounds, availability of the complete 2^6 matrix element in Monte Carlo is mandatory, see eg. discussion in |35 1. The 
total cross-sections are given in Table|7|for different definitions of oqed, o-qcd couplings. 

As an example, Fig.|2|illustrates spin correlation effects in the top-pair production and decays, namely asymmetry in the 
correlations between lepton and antylepton direction in the rest frame of top-quark, for events generated with 2 — > 6 matrix 
element. Such correlation is absent if only 2^2 matrix element is used for events generation, followed by the independent 
decays of each top-quark. 
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Table 7: AcerMC production cross-sections for the gg, qq {WWbb -^) f f f fhh process. 
The 14 TeV centre-of-mass energy and CTEQ5L parton density functions were used 
along with different definitions for the aqED, oiqcd (as in native AcerMC, default 
PYTHIA 6.2 and HERWIG 6.51 



Factorisation scale 


native AcerMC 


CtQED, OIQCd(IL) 

asin PYTHIA 6.2 


CCQBD, OIQCD (2i) 

as in HERWIG 6.5 




gg {WWbb -^)ffffbb 


Qqcd = (2 • m") 


400. [pb] 


390. [pb] 


330. [pb] 


Qqcd = E (pi," + m?)/4 


450. [pb] 


450. [pb] 


380. [pb] 


Qqcd = E (Pt')/2 


550. [pb] 


550. [pb] 


460. [pb] 


Qqcd — * 


355. [pb] 


350. Ipbl 


300. [pb] 




qq ^ (WWbb -^)ffffbb 


Qqcd = (2 • m^) 


63. [pbj 


62. Ipb] 


53. Ipbl 


Qqcd = E ij>V 


69. [pb] 


69. [pb] 


58. [pb] 


Q%CD = E (pi.')/2 


78. [pb] 


78. [pb] 


65. [pb] 


Qqcd — * 


59. [pb] 


58. [pb] 


50. [pb] 



Table 8: AcerMC production cross-sections for the gg,qq — > (tt —^)ffbffb). The 14 TeV 
centre-of-mass energy and CTEQ5L parton density functions were used along with 
different definitions for the uqed, olqcd (as in native AcerMC, default PYTHIA 6.2 
and HERWIG 6.5;. 



Factorisation scale 


OIQED, ocqcdC^L) 
native AcerMC 


aQED, aQCD(iL) 

as in PYTHIA 6.2 


aQED, olqcdC^L) 
as in HERWIG 6.5 




gg ^ (tt -*)ffbffb 


QlcD = (2mt) 


370. Ipb] 


365. [pb] 


310. [pb] 


QlcD = E (Pt''+">?)/2 


425. Ipb] 


420. Ipb] 


355. [pb] 


QlcD = E {pV)/2 


512. Ipb] 


510. Ipbl 


425. [pb] 


Qqcd = « 


330. [pb] 


320. Ipbl 


280. Ipb] 




m - It -)fl!'lfi> 


Qqcd = (2 • m") 


63. Ipbl 


62. [pb] 


53. [pb| 


QlcD = E {pV + 


69. [pb| 


68. [pb] 


58. [pb| 


QlcD = E (Pt')/2 


78. [pb| 


78. [pb] 


65. [pb| 


Qqcd — ^ 


59. [pb] 


57. [pb] 


50. [pb] 
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Fig. 2: The correlations between cos (azimutal angle) of lepton and antylepton from tt — > 
iDbih'b decays measured in the rest frame of the top-quark with respect to the anty- 
top quark direction. Left plot is for gg {WWbb —>-)ffffbb process, right plot for 
qq (yVWbb —>^)ffffbb process. 



4.8 The gg,qq tttt process 

This process, in spite of having a very small cross-section at LHC energies, contributes as reducible background to the 
HH WW WW production at low masses, 1131 . Availability of the complete Monte Carlo generator is mandatory to give 
reliable predictions of theis background and to optimise selection criteria. 



Table 9: AcerMC production cross-sections for the gg, qq — > titi process. The 14 TeV centre- 
of-mass energy and CTEQ5L parton density functions were used along with different 
definitions for the agED, olqcd (as in native AcerMC, default PYTHIA 6.2 and HER- 
WIG 6.5j. The niH = 120 GeV and mt ~ 175 GeV were used for calculating the 
Qqcd- 



Factorisation scale 


OiQED, OLQCDi^L) 

native AcerMC 


CtQED, aQCD(li) 

as in PYTHIA 6.2 


OiQED, aQCD{2L) 

as inHERWIG 6.5 




gg — > tttt 


Qqcd — ^ 


26.5 [fb] 


24.4 [fb] 


19.3 [fb] 


Q'ocD = E (Pr' + ™')/4 


75.7 [fb] 


73.8 [fb] 


53.4 [fb] 


Q%CD = Z (Pt')/4 


94.7 [fb] 


93.2 [fb] 


66.2 [fb] 


Qqcd = C™* + mff /2)^ 


89.5 [fb] 


87.8 [fb] 


62.9 [fb] 




qq — > tttt 


Qqcd = ^ 


0.5 [fb] 


0.5 [fb] 


0.4 [fb] 


Qqcd = Y. {pi-'" + mf}/i 


1.2 [fb] 


1.2 [fb] 


0.9 [fb] 


QlcD = E iPT'")/4 


1.5 [fb] 


1.5 [fb] 


1.0 [fb] 


Qqcd = ("'■t + ™ff/2)^ 


1.5 [fb] 


1.4 [fb] 


1.0 [fb] 
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4.9 The control channel processes 

The set of simple 2-^2 control channel processes was added to AcerMC in order to provide a means of consistency and 
cross-check studies. Although these processes are already implemented in PYTHIA and/or HERWIG (except the gg — > 
WbWb one), the availability of the native implementations is supposed to offer a more consistent control of generation 
parameters when performing e.g. the comparison of parton shower PYTHIA/HERWIG produced additional pair of heavy 
quarks with the exact leading-order matrix elements implemented in the core group of the AcerMC 2^4 processes. 

To the list of control channel processes we have added also 2 — » 4 process, the gg WbWb process, as a control channel 
for the 2^6 process gg ffbffb. The 2^4 process we consider as very usefull for studying in more detail the resonant 
and complete WbWb production at LHC. 

The benchmark results, given in TablellOl are obtained for simulated pp collisions at 14 TeV centre-of-mass energy. The 
total cross-sections are listed for different definitions of uqed, ctqcD couplings and different definitions of the energy scale 
Q\ 

Table 10: AcerMC production cross-sections for the qq — > Z/^* —> U, gg,qq — > ti and 
qq ^ W ^ Iv with single flavour Z/j* i£ and/or W decays. The 14 TeV 

centre-of-mass energy, CTEQ5L parton density functions and different definitions for 
the ctQED, OiQCD (os in native AcerMC, default PYTHIA 6.2 and HERWIG 6.5) were 
used. The threshold mu > 60 GeV was used in the event generation. In PYTHIA 
6.2 the setting MSTU(1 15)=2 was used to set the lower limit in oqcd evolution 
to 4 GcV^ as done in the native AcerMC implementation. 



Factorisation scale 


native AceiMC 


OiQED, aQCoilL) 

as in PYTHIA 6.2 


aQED, oiqcdC^L) 

as in HERWIG 6.5 




qq ^ ^/7» ^ U 


= s 


1620 [pb] 


1630 [pb] 


1630 [pb] 


= (E (pt') + m|)/2 


1550 [pb] 


1560 [pb] 


1560 [pb] 


= E (Pr')/2 


1260 [pb] 


1260 [pb] 


1260 [pb] 




1630 [pb] 


1630 [pb] 


1640 [pb] 




gg tt 


= s 


365 [pb] 


360 [pb] 


310 [pb] 


= i:(pv' + ™?)/2 


430 [pb] 


420 [pb] 


355 [pb] 


= ■£ (pV')/2 


595 [pb] 


590 [pb] 


490 [pb] 


= {2mtY 


320 [pb] 


315 [pb] 


270 [pb] 




qq — > tt 


= s 


bl. [pb] 


61. [pb] 


52. [pb] 


= Y. (pV" +™?)/2 


69. [pb] 


68. [pb] 


58. [pb] 


= E (pt')/2 


86. [pb] 


85. [pb] 


71. [pb] 


= (2mt)^ 


57. [pb] 


56. [pb] 


48. [pb] 




qq~,W 


= s 


17200 [pb] 


17230 [pb] 


17310 [pb] 


Q"- = (E (Pt') + a4')/2 


16480 [pb] 


16490 [pb] 


16460 [pb] 




12920 [pb] 


12920 [pb] 


13020 [pb] 


Q"- = Ml- 


17360 [pb] 


17380 [pb] 


17300 [pb] 




gg WbWb 


= S 


370 [pb] 


365 [pb] 


310 [pb] 


Q' = E (Pr' +™?)/2 


430 [pb] 


425 [pb] 


355 [pb] 


= E (pt')/2 


525 [pb] 


520 [pb] 


435 [pb] 




330 [pb] 


320 [pb] 


275 [pb] 
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5 Monte Carlo algorithm 

The conceptual motivation leading to the present implementation of AcerMC was to exploit the possibility of dedicated 
matrix-element-based generation interfaced to a more general event generator, called supervising event generator, which is 
subsequently used to complete the event generation procedure. 

The goal of the dedicated matrix-element-based part is to efficiently generate complicated event topologies using native 
(multi-channel based) phase space generation procedures. The strategy is based on the understanding that a case-by-case 
optimisation is in complex cases of phase space topologies preferable to an universal algorithm. Given that phase-space is 
optimised on a case-by-case basis, an user-defined pre-selection for the generated regions of the phase-space is not imple- 
mented. Due to the fact that the 2 ^ 4 and 2 ^ 6 matrix elements, provided by the MADGRAPH/HELAS |8| package, 
contain full massive treatment of the final state particles, there are no explicit divergences present for implemented processes 
and AcerMC can indeed cover the full (kinematically allowed) phase space of the processes at hand. 

The matrix-element-based part uses aQcoiQ^) and aQEoiQ^) couplings and mass spectra, as calculated by the super- 
vising event generator, to insure the full internal consistency in treatment of the event itself. Optionally, the native oiqcd (Q^) 
and oiqed{Q^) definitions can also be invoked. 

The generation chain is built from the following steps: 

• The PYTHIA 6.2 or HERWIG 6.5 interfaces to the library of the structure functions PDFLIB 8.04 |4T1 are used to 
calculate convolution of the partonic density. 

• AcerMC modules produce unweighted hard-process events with colour flow information and pass them to the super- 
vising generator PYTHIA 6.2 or HERWIG 6.5 as an external event. 

• The generated events are then further treated within PYTHIA 6.2 or HERWIG 6.5 event generators, where the frag- 
mentation and hadronisation procedures, as well as the initial and final state radiation are added and final unweighted 
events are produced. 

The AcerMC efficiency* for generating unweighted events, using the implementation of the phase-space generation 
discussed below, is summarised in Table fTTI A certain (very small) fraction of events is further rejected in the shower- 
ing/fragmentation procedures of the supervising generators. 

In the following we will briefly describe the key points of the implemented AcerMC modules and developed algorithms: 
matrix element calculations, n-fermion phase-space generation based on the modified Kajantie-Byckling methods| 54 1, the is- 
sue of the s-dependent width and mass threshold effects for resonances and finally, the modification of the VEGAS algorithm. 

5.1 The Matrix Element Calculation 

The squared matrix elements of the processes were obtained by using the MADGRAPH/HELAS jsj package. They take 
properly into account the masses and helicity contributions of final states particles, incoming quarks are considered as mass- 
less. The particle masses, charges and coupling values that are passed to the code derived with the MADGRAPH package 
are calculated from functions consistent with the ones used in supervising generators (PYTHIA/HERWIG). This allows to 
preserve the internal consistency of the event generation procedure. In particular, the (constant) coupling values of and 
oqed were replaced with the appropriate running functions that were either taken from the interfaced generators or provided 
by the AcerMC code according to user settings. Slightly modified MADGRAPH/HELAS allowed for obtaining colour flow 
information of the implemented processes. 

The sets of the MADGRAPH/HELAS coded diagrams, for each of the implemented processes, are collected in Appendix 

A. 



*Note that efficiency is energy scale dependent and phase-space optimisation is done individually for each choice. So it might vary for the same process 
but different choices of the energy scale definition. 
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Table 11: Efficiency for the generation of unweighted events with the default definition of the 
energy scale, ACES2=1 (see Section Wy\f or details). For generation of the qq, gg — » 

Z/"f*{^ i£)bb and qq,gg Z/^*(-^ i£)tt events threshold mil > 60 GeV has 
been used. The f = e, ii,t, q,b. 



Process 


Description 


Internal AcerMC efficiency 








[1] 


gg — > ttbb 


36.3 % 


[2] 


qq ttbb 


29.7 % 


[3] 


qq^W{-^ vt)bb 


35.2 % 


[4] 


qq W[-^ vt)tt 


30.0 % 


[5] 


gg Z/-i*{-^ U:)bb 


42.8 % 


[6] 


qq^ Z/"f*{-^ ££)bb 


35.1 % 


[7] 


gg^Z/Yi-^ ff,i^v)ti 


47.0 % 


[8] 


qq^Z/Yi^ ff,y^)ti 


42.6 % 


[9] 


gg ^ [Z/Wh* -^)ttbb 


9.3 % 


[10] 


qq {Z/Wh* -^)ttbb 


32.4 % 


[15] 


gg tttt 


48.0 % 


[16] 


qq tttt 


50.2 % 




2^6 




[11] 


99^{tt^)ffbffb 


14.2 % 


[12] 


qq^ {tt^)ffbffb 


12.0 % 


[13] 


gg ^ {WWbb -^)ffffbb 


18.2 % 


[14] 


qq ^ (WWbb -^)ffbffbb 


4.2 % 




Control processes 




[91] 


qq Z/j* ££ 


68.4 % 


[92] 


99 tt 


65.6 % 


[93] 


qq — > tt 


62.1 % 


[94] 


qq^W -*vl 


69.4 % 


[95] 


gg (tt ^)WbWb 


40.2 % 



5.2 The Phase Space Generation Procedure 

The general objective in simulation of physics processes for the LHC environment is to improve the integration of the differ- 
ential cross-section using Monte-Carlo sampling methods'. The sampling method used should aim to minimise the variance 
of the integral as well as maximise the sampling efficiency given a certain number of iterations and the construction of the 
sampling method itself should aim to be sufficiently general and/or modular to be applicable to a wide range of processes. 
Writing down a (process) cross-section integral for LHC type (hadron-hadron) collisions: 

cr = y faixi, Q'^)fb(x2, Q^) (-2j)t!-4(-2g) '^^^ 

where fa,b(x, Q^) represent the gluon or (anti)quark parton density functions, |A^n|^ the squared n-particle matrix element 
divided by the flux factor [(27r)^"^*2s] and d^n denotes the n-particle phase space differential. The quantity s = xi X2 s is 
the effective centre-of-mass energy, and the sum ^ runs in case of quark-antiquark incident partons over all possible quark- 
antiquark combinations (a, b = u, d, s, c, u, d, s, c). In case of gg initial state the sum has only one term with a = b = g. 

'For a nice discussion on the topic see e.g. I49II50I . . . 
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It is often convenient to re- write the differential cross-section in the form: 

a.b 



(27r)3"-4(2s2 



■ dy ds d^„ 



(2) 



with the new (rapidity) variable given by y = 0.5 log(xi/x2). The n-body phase-space differential d$n and its integral 
$n depend only on s and particle masses rrii due to Lorentz invariance: 



<E>„(s,mi,m2, 



d$„(s,mi,m2, . . . ,m„) = / 5* (pa + Pi) - ^ p,: d''pi(5(Pi - )e(pi ), (3) 



i=l / 1=1 



with a and b denoting the incident particles and i running over all outgoing particles i = 1, . . . ,n. What one would like to do 
is to split the n-body phase parameterised by 3n-4 essential (i.e. non-trivial) independent variables into manageable subsets 
(modules) to be handled by techniques which reduce the variance of the result and/or the sampling efficiency (e.g. importance 
sampling|40 1 or adaptive integration like VEGAS jQj or FOAM 1441 ). Stating this in formal terms, the above Equation|2|should 
be transformed into an expression like: 



/ 4 \ 

n i. 

n / 



\ 



I 



I 



J=l _ 



t \ 



V 



J-a\ \ ^ a;i/a(xi, Q^) a;2/i,(x2, Q^) 



|Mn| 



(27r)3"-*(2s2) 



dy ds (4) 



where one integrates over Mandelstam type (Lorentz invariant) momenta transfers Si , tj and space angles Qk = (cos ■&k,<j)k) 
within the kinematically allowed limits (3n-4 variables in total) with the term IJn j denoting the Jacobian of the transformation. 
If one would then decide to introduce importance sampling functions in order to reduce the peaking behavior of the integrand 
1401 . the integrals would take the form: 

s+ s + 



dsi 



gi(si) 



dsi, 



(5) 



where the importance sampling function gi is probability density function normalised in the integration region [s; , 



i 



gi{si)dsi 



(6) 



which exhibits a similar peaking behavior as the integrand. Formally, one then inserts the identity: 







gi{si)dsi 



\ 



dn 



I 



into the integral and then derives the unitary sampling prescription: 



1 / 
dn I 5 



gi{si)dsi 



\ 



1 "i 1 
= [ dr^ [ 5{si- G"^(rO) —j-^dsi = / J'^' 

g^{si) J J ^ gi{Si) J gi{G ^CO) 





(V) 



(8) 



which formally means that the Si values are sampled from the interval according to the gi (si) distribution by using (pseudo-)random 
variable li together with the gi(si) cumulant G(si) — ^"l gi{si)dsi with its inverse G^^. The unitarity of the algorithm states 

that each trial ( ri value) produces a result (i.e. a corresponding si value distributed according to gi(si)). 
Performing such substitutions on all integration parameters would give as the cross-section expression; 



3n-4 j. 

CT = ]^ dri 

j— 1 r\ 



f{ri,r2...) 
g{ri,r2 ■ ■ .) 



(9) 
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where the integrand would (hopefully) have as low variation as possible at least for a subset of contributing Feynman dia- 
grams** . To improve the sampling method further, the ri (pseudo-)random variables can be sampled from adaptive algorithms 
of the VEGAStype O. 



9 b a b Pb Pi 




Fig. 3: A representative Feynman diagram describing a 2 6 process gg ^ tt ^ bhW^W — > 
bblueq\q2 and its decomposition into a set of 2 2 t-channel and s-channel sub-processes. 

A representative Feynman diagram describing a 2 ^ 6 process is shown in Figure|3| As one can see, the process can be 
split in several consecutive branchings, this approximation is often used in matrix element (probability amplitude) calculations. 
It seems rather obvious that any Feynman diagram can be split in a series of horizontal and vertical branchings that one can 
denote as s-type and t-type(u-type) using the analogy with the Mandelstam variables. What one would like to do is thus to 
modularise the phase space in the form of sequential s- and t-type splits. 

The s-splitting of phase space is relatively easy to do and has as such been used in many instances of Monte-Carlo 
generation (e.g. FermiSV |37|, Excalibur |38|, Tauola |5| etc.); the t-type branchings (often tagged as multi-(peri)pheral 
topologies) have in contrast generally been calculated only for specific cases (e.g. for 3 or 4 particles in the final state I37II38I ). 
As it turns out, the problem of several massive particles in the final state has already appeared more than 30 years ago when 
several hadrons (e.g. pions) have been produced in (comparatively low energy) nuclear interactions. At that time Kajantie 
and Byckling 1 54!j have derived the formulae for simulating any sequence of s- and t- type branchings which, with some 
modifications, can also be applied to the EW and QCD processes involving heavy quarks and/or massive bosons at LHC. 



^The 'modularisation' can be performed for several topologies at the same time and multi-channel techniques can be applied. 



B. P. Kersevan,E. Richter-Was: 



AcerMC 2.0 Event Generator 



22 



5.3 Modified Kajantie-Byckling Formalism 
5.3.1 The s-type Branching Algorithms 




Fig. 4: A diagram of a generic 2 — > 2 s-channel process. 

The s-splits are the simplest method in the KB formahsm. For the sake of completeness one should start with the definition 
of the two-body phase space integral (c.f. Fig|3J: 



$2(5, mi, 7712) = J d'*'pid'^p25{pi - mi)5{p2 - m2)5*{p - Pi - p2)0{pi)0{p2), 



(10) 



with the incoming momentum sum p — (pa + Pb),p^ ~ s and outgoing momenta pi,2,Pi,2 = nii,2- The phase space 
integral is Lorentz invariant (as one can observe in the above Equation where it is written in a manifestly Lorentz invariant 
form). Subsequently, due to Lorentz invariance, the integral is necessarily a function of the Lorentz scalars s, mi and m2 only. 
The step function product O(pi)0(p2) is the explicit requirement of the positiveness of the energy terms in pi,2 while the 
delta functions represent the on-shell conditions on p 1,2 and the total momentum conservation. 

The integral can be transformed into a more compact form by integrating out the spurious variables; one thus first integrates 
over d''p2 and chooses the centre-of-mass system (CMS) as the integration system of reference with p — {y/s, 0, 0, 0) and 
then evaluates the integrals over p? and : 

$2(3,7711,7712) = J d^pi5{pl - ml)5{{p - pif ~ ml)0{p1) (11) 
d'pt 



^S{s + mt- 2^sEl - : 



s + nil ^ m,2 



2El 

^ / P^dEld^S (eI - 2^ 

Pl(s, 7771, 7772 



4^ 



dfii, 



with the stars explicitly denoting the values in the centre-of mass system. The first integration simply sets pi = \/ (Pi)^ + nij - 
and the second integral leads to the well known relations for the energy: 

and momenta sizes: 



Pi = bi! = -^^ ^^/l , P2=Pi (13) 

of two particle production. The A(s, mi, m|) denotes the Lorentz invariant function: 

A(s, 777?, 7772) = (s - (7771 + 7772)^)(s - (777i - 7772)^) (14) 

and thus explicitly contains the phase space cutoff, i.e. the requirement that the available CMS energy ^/% should be bigger 
than the mass sum ^/% > (mi + m2). Note that the integration was so far done only over the spurious parameters, leaving the 
polar and azimuthal angle of the pi particle as the two independent parameters dfl* = dcosO*d(p* . The integral becomes 
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trivial to sample in case the outgoing particles can be approximated as massless (the 'boost' factor lambda transforms to 
unity). As already claimed, the latter approximation is however often unjustified when studying processes representative for 
the LHC environment. 

Kajantie and Byckling f54\ introduced the recursion and splitting relations for the n-particle phase space $n(s) given by 
Eq.|3] The recursion relation can be derived by defining the momentum sum: 



J = l 



Subsequently one can interpret p = kn and s — from Eq.|3| One continues by introducing the identities: 

1 = j dM^.sikl^, - M^i)e(fcn-i) 



and 



1= / d'^kn~\5'^{p — k 



n-l — Pn 



(15) 

(16) 
(17) 



(18) 



into the integral of Equation|3l separating out the arguments containing kn-i andp„ terms one obtains: 

<l>„(M^,mi,m2, . . . ,m„) = J dM^_i x 

X |y d'^kn-id'^pn5{ki_i - Ml^_i)S{p'i - mi)5'^{p - kn-i - pn)Q{k°_i)Q{p'l)^ x 

X $„_i(M^_i,mi,m2, . . . ,m„_i), 

where the remaining pi terms form the (n-l)-particle phase space integral $n-i (M^-i, mi , m2, . . . , mn-i) and the terms in 
curly brackets give a two particle phase space term (c.f. Eq. I12> : 



$„(M^,mi,m2, . . . ,m„) = J dM^_i$2(M^, M„_i, m„)$„_i (Mi5_i , mi, m2, . . . , mn_i) 



(19) 



,mi,m2, . . . ,mn~i) 



dMl_^ 



8M2 



, mi, 7112, • . • , m„_i ) 



with the integration limits on Mn_i following from its definition in Eg. 1151 It has to be emphasized that the angles in dO.* are 
each time calculated in the centre-of-mass system of ki with the invariant mass Mi. The resulting recursion relation is clearly 
of advantage when describing cascade decays of particles fe„ fc„_ip„ ^ kn~2,Pn,Pn-i ^ . . .; it also proves that the 
n-particle phase space of Eq.|3|can be reduced into a sequence of two-particle phase space terms, as shown in Figure|5| 



Pn Pn-1 Pn-2 



Pi+1 Pi 




Fig. 5: The diagrammatic representation of consecutive s-splits. 
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It can further prove of advantage to loosen up the splitting terms of Eqns. ll6ll7l so that instead of summing to n-1 one 
groups an arbitrary set of £ particles: 

1 = J dM?5{kf - Mf )e(fc?), (20) 

/n 
d*kiS\p-ki- J2 P^)' (21) 

which, when repeating the procedure in recursion relation of Ea. fT9l results in an expression: 

(A--fi + i-m,^l)^ 

$„(M^,mi,m2, . . . ,m„) = J diWf $„_;+i(M^, M;, mj+i, . . . , m„)<E>i (Aff , mi, m2, . . . , m;), (22) 

and thus effectively splits the phase space into two subsets, equivalent to introducing an intermediate(virtual) particle with 
momentum ki . 

The number of splitting relations and the number of particles in each group as given in Eq. 12 II can be chosen in any 
possible sequence, thus meaning that the grouping sequence is arbitrary and can be adjusted to fit the topology in question.' 

At this point some modifications were introduced to the algorithm in order to adapt it to the specifics of the processes 
expected at the LHC. Kajantie and Byckling namely assumed that the generation sequence would be 'down' the cascade (i.e. 
by sampling first a Mn value, then M„_i value etc. . . as is indeed most often done in Monte-Carlo Generators). This might 
however not be optimal in the LHC environment since the available centre-of-mass energy for the hard process (s) can vary 
in a wide range of values (c.f. Equation |2j and has to be sampled from a distribution itself. The shape of the distribution 
function for s is expected to behave as a convolution of the peaking behavior of all participating invariant masses times the 
parton density functions (c.f. Eq. |2j; it subsequently seems to be more natural (and efficient) first to sample the individual 
propagator peaks and then their subsequent convolutions. Furthermore, by generating the invariant masses 'up' the cascade 
(i.e. first M2, M3 . . . Mn and finally s) the kinematic limits on the branchings occur in a more efficient way (bound on the \/A 
values, see Eauations ll4l and l24l . which is very convenient since in the LHC environment no stringent generation cuts should 
be made on the inherently non-measurable s as it cannot be accounted for by an analogous cut in a physics analysis. 

A necessary modification of the algorithm would thus be to reverse the generation steps by starting with the last pair(s) of 
particles. In terms of integration (i.e. sampling) limits this translates into changing the limits of Eg. 1191 

$„(M^,mi,m2, . . . ,m„) = (23) 

Jx{Ml,Ml_,,ml) r 

dML.^ — — y^"" 

(Af„_i-m„_i)2 



X ... J dM!^ Jdn: 



{mi+m2)^ 

which accommodates the mass generation sequence: kn — > + pn ■ ■ ■ (i.e. first sample M^__i, then A'In-2 etc. . . ), 



Suggestions of 1 54 1 on how to pick random number sequences will not be used since one might like to couple this method with an adaptive algorithm to 
improve the sampling efficiencies. 
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into 



{M„-m„)2 

$„(M^,mi,m2,...,m„) = J dMl^^^ — ^^^75— j dQ,l (24) 

(M„_2+™.„-l)^ 



"8M| 



X / ^2 / '^"--1 

(M„_3 + m„_2)2 

X ... y dMi^v — _ I 

(JVIi_l+mi)2 
(A-f„-E"=3"ij)^ 



^ 8iw| 7""^' 

(mi+m2)-^ 

where one first samples the mass M2, M3. . . M„_i in the appropriate limits. 

In some topologies symmetric cases of mass generation can appear (as shown in Figure|3} where the integration sequence 
is ambivalent (e.g. in Figure|3|the ambivalence is which top quark invariant mass to generate first. . . ) and after a choice 
is made (since one of the two cases in the symmetric pair has to take precedence) the procedure itself remains not entirely 
symmetric. Detailed studies have shown that it proves useful to include all permutations of such ambiguous sequences into 
the MC algorithm in order to 'symmetrise' the solution and thus make it easier to process by further additions (e.g. adaptive 
algorithms). 

5.3.2 The t-type Branching Algorithms 




Fig. 6; A diagram of a generic 2^2 t-channel process. 



The t-splits are a specialty of the KB formalism due to the advanced calculation of the limits on the (massive) t- variable. 
The formalism can be introduced by observing that in case of a pa + Pb — > Pi + P2 scattering the momentum transfer is 
characterised by the (Mandelstam) variable t — (pi — pa)'^ (c.f. Fig|SJ. It is thus sensible to replace the = dcos 61dip\ 
integration in the two body phase space integral of Eq. E|with integration over the t variable. Writing the definition of t in 
the centre-of-mass system one gets: 

t = q^ = {Pa^Pif (25) 

2 



mi - 2ElEl + 2p*pi cos 61 



and hence: 



dt^2plpld cos e* (26) 
Using the latter substitution together with Eq. 1121131 and the analogue for pa: 

\/ X(s, m?,,m?) 

Pa = ^ (27) 
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one obtains in place of Eq. 1121 



$2(s,mi,m2) = 



Pi(s,mi,m2) 



4^ 



J dnl 



(28) 



dt dip* 



4,^\{s,ml,ml) 



+ 2TT 

dt j dtp* 





With the integration variable change the integration domain changes from [—1, 1] fordcosS* to [t ,t"'"] for the dt integration. 
The limits are obtained by inserting the cos 9* limits into EQuation l26l 

t^^ml+ml- 2E:eI ± 2p:pt (29) 

or in the Lorentz invariant form (c.f. Eq. 1121131 : 



,± _ 2,2 (s + m„ -mft)(s + mi -ma) 

I — "la + '"l n 



(30) 



^\{s,ml,ml)\{s,ml,ml) 
2s 



As a step towards generalisation one has to note that the kinematic limits t can also be derived from the basic four- 



particle kinematic function Gfa.v.z.t<.v.w) l55ll54l . where the function G can be expressed as a Cayley determinant: 

G{x,y,z,u,v,w) = - 
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The kinematic limits on t are in this case given by the condition 



G{s,t,m2,ma,mt,mi) < 0, 



(31) 



(32) 



it should be noted that the above condition gives either limits given a fixed value of s or equivalently limits given a 
fixed t value. 

In search of a recursion relation involving t-variables one can note that in Eq. I19l the angle in cos is equivalent to the 
scattering angle in the centre-of-mass system of the reaction Pa + Pb ^ kn-i + Pn and thus given by: 



tn-l = {pa — fcn-l) 

= ml + M^_i - 2E*k'},*_i + 2p*fc*_i cos0*„i 



(33) 



with the tj7_j limits expressed by: 

G{Mlt„^i,mi,ml,mlM^_i) <0, (34) 

and the Pa given by Eq. 1271 In order to produce a more general picture it can further be deduced that the next angle in the 
recursion is the scattering angle of the subsequent process pa + (pb — Pn) ^ kn-2 + Pn-i in the centre-of-mass 

system of fc„_i; the (pb — Pn) = qn-i is in this case considered as a virtual incoming particle with momentum qn-i (c.f. 
Figure 13 • 

It immediately follows that for a general process pa + Qi+i —> ki + Pi+i with: 



Pb- ^ Pj ^ Pa - ki 

3 = i + l 



11 2 
Qi — ^i] Qn — in = TTii^ 



a general expression for ti becomes in the centre-of-mass frame of ki+i: 

— {pa ~ ^i) 



(35) 



(36) 



ml + M!- 2i5:('+i)fcf (^+^) + 2p*J'+'h*'-^+'^ cose* 
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Fig. 7: The diagrammatic representation of the method apphed in translating the multi-(peri)pheral splits 
into a 2 — > 2 t-channel configuration. 



where momenta in centre-of-mass frame of ki+i, denoted with the superscript *(i + 1), are given by: 



.(i + l) 



2M,+i 



(37) 



(38) 



and the corresponding energies fc^ * and Ea can simply be obtained by using the analogues of Equations ! 121 IBI or the 
usual Einstein mass-energy relations directly. The corresponding if limits given by: 

G{Mf+i,ti,mj+uml,U+i,Mf) < 0, (39) 

and the recursion relation of Ea. ll9l becomes: 

$„(M^,mi,m2, . . . ,m„) = 

{M„-m„) 

: J dip*n j dtn-1 "I>n-1 ,mi,m2, . . . ,m„_i) 



(40) 



As already argued the resulting set of (si = M?, ti) can again be sampled in any direction with respect to the cascade 
by applying the appropriate change in the integration limits (c.f. Eq. 1191 and 1241 . The recommended approach (i.e. the 
introduced modification of the algorithm) is again to first sample the invariant masses in the reverse cascade direction (i.e. in 
the sequence M2, M3, . . . , Mn) and then the ti values within the limits calculated from Eg. l34l down the cascade (i.e. in the 
order of U-i, tn-2, • • • , ti). 

To sum up, it has been shown that using the Kajantie-Byckling formalism the phase space for any topology can be split 
in a set of s-type and t-type 2^2 branching steps (modules) given by recursive formulae of EQuations l24l and l40l 
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5.4 Propagator Sampling 

A well known theoretical issue is that one can expect the most prominent peaks in the differential cross-section of a specific 
process in the phase space regions of high propagator values in the corresponding probability density. Consequently, in the 
scope of complementing the modular structure of the derived Kajantie-Byckling based phase space sampling, new approaches 
were also developed concerning the numerical sampling methods of the relevant kinematic quantities. 

In order to get small variance in the Monte Carlo procedure one would thus like to include the appropriate peaking 
dependence of the relevant momentum transfers in the importance sampling function. It however turns out that since the 
momenta transfers q participate also in the propagator numerators (typically in p^q^/q^) and since in process of interest 
one mostly finds several Feynman diagrams contributing to the final probability density, thus causing interferences, it is very 
difficult or even impossible to estimate the exact power of momenta transfers in the sampling functions for different propagator 
peaks. In other words, the probability density dependence on the momentum transfer q^ can in general be approximated with 
the dependence l/(q^)'^ where the best value of v must be determined separately (on a process by process basis). 

In view of the latter, general formulae have been developed for sampling the x^^ shape 1.37. .56J : Given a pseudo random 
number r G [0, 1] and limits x € [x-,x+\ the value x distributed as x~" is obtained from the formulae in Eq.|8|as: 

X = [x:''+'-(l-r)+xI"+'-r]"^; u ^ I; (41) 
X = v = \. (42) 



Using the analogous (unitary) approach a recipe for resonant (Breit-Wigner) propagator contributions of the type: 

with s £ [s_ , s+] and a pseudo random number r G [0, 1] is available by the prescription: 

s = + Mr • tan[(u+ - u_) ■ r- + tt_] (44) 
atan — — (45) 



± 



Following similar arguments as for the non-resonant propagators one can surmise that the best sampling function for 
resonant propagators could in general be a Breit-Wigner shape modified by a factor ^ ,v ^ [0, 1]. As shown in the following 
section it was found that a shape: 

works quite well for a set of processes and a corresponding sampling recipe was developed. In addition, studies in f57"| show 
that a resonant -^sx Breit-Wigner shape: 

should be expected in a range of decay processes. Detailed studies have shown that it is in general better to introduce a 
s",!^ G [0, 1] dependence even if it over-compensates the high mass tails of the corresponding differential cross-section 
distribution since this provides an overall reduction of the maximal weight fluctuations in the Monte-Carlo event generation 
procedure. 

5.4.1 Breit-Wigner Function witli s-dependent Width 

In some topologies of the processes involving or bosons, a bias of the matrix element towards large values in the 
high region is evident, which in turn means that a more accurate description of the tails of s^y^ distribution is needed. 
Consequently, the Breit-Wigner sampling function was replaced by'": 



To our knowledge this implementation is original and done for the first time in AcerMC. 
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which is proportional to the (more accurate) Breit-Wigner function with an dependent width (W in the above formula 
denotes either a or a Z'^ boson). 

In order to implement a unitary algorithm (an algorithm that produces a result for every trial, i.e. there is no rejection) 
of value sampling on the above function one first has to calculate the normalisation integral (cumulant) and then its inverse 
function. Introducing a new variable rj = (sj^,^ — A/^)/ {Mw ^w) the integral of the above function can be expressed as: 

where the upper integral limit is left as a free parameter. The integral thus gives a function: 

F{V) = { ■ atan(r,)| + | i ■ log(^^ + 1) | , = Fi (,7) + F^iv) (50) 

with -F(?;max) — F{riinin) defining the normalisation. One of the undesirable features is that the function F{rj) does not have 
a (simple) analytical inverse, which is a prerequisite for unitary sampling. Taking a closer look at the two above expressions 
one can quickly spot another undesirable feature, namely that the second term in the Equation 1491 is an odd function of rj, 
which after the integration gives an even term -^2(^7) in rj in Eguation lSOl In other words the second term alone is neither a 
non-negative function nor does it have an unique inverse - one has to deal with a negative probability. A reasonably elegant 
solution to this problem has been developed and implemented here: 

• One samples values of 77 by using only the first term of the above expressions (the usual Breit-Wigner function). 

• One then re-samples the obtained value of r) using the full expression of Eauation l49l If 77 is less than zero the value is 
mapped to —77 with the probability given by Eauation l49l 




Fig. 8: Left Comparisons of the two functional terms of Eq. [49\ to JiW s{ti) given by Equation 
[48\ Note that the scaling factor A is chosen in view of making the contributions more 
transparent; it is much too small compared to the real case ofW^/Z^ bosons. 
Right Comparison of the (normalised) distributions of differential cross-section for the 
process qq Wbb (dashed) and sampling functions (solid line) with respect to the 
variables obtained by importance sampling, as described in the text. 



Why this works can quickly be deduced by looking at the Figure|8| At negative values of r) the second term of Equation 
1491 gives a negative probability in the region 77 < 0, i.e. using a simple Breit-Wigner (Cauchy) probability function too many 
events are generated in this region. Correspondingly, since the second term of Eq. I49l is an odd function, exactly the same 
fraction (distribution) of events is missing in the region 77 > 0. By mapping events with 77 < over the r) = G axis one thus 
solves both problems at the same time. Using the above re-sampling procedure the whole approach remains unitary, i.e. no 
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events are rejected when there are no hmits set on the value of rj or they are symmetric |r7min| = »7max- In the contrary case, a 
small fraction of sampling values is rejected. 

After some calculation the whole unitary procedure can thus be listed as follows: 

• Calculate the kinematic limits 77miii and 77max. 

• Calculate the normalisation factors Ai — fi(77max) — fi(»7min), A2 = F2(r7niax) — f2(r7min) and As = Ai + A2; 
the term A2 can actually be negative and thus does not represent proper normalisation. 

• Obtain a (pseudo-)random number pi . 

• If pi < A2/AS then: 

- Obtain a (pseudo-)random number P2; 

- Construct rj as: 

X = A2 ■ P2 +F2(77n,in), 



= Vie 



which is the inverse of the (normalised) cumulant {F2{rf) — -F2(?;min))/A2. 

- Note that the condition pi < A2/AS can be fulfilled only if A2 > 0, which means that 77max is positive and 
greater than 77min. 

Conversely, ifpi>A2/As then: 

- Obtain a (pseudo-)random number P2; 

- Construct 77 as: 

X = Al ■ p2 +Fl(r?niin), 

,Mw^w „N 
^ = tan(— ^-X) 

which is the inverse of the (normalised) cumulant {F\(ri) — Fi(r;min))/Ai. 

- If the obtained 77 is less than zero then calculate the normalised probability densities: 

P = ^ r _J_ 

^ Ai ^ MwTw l + v'^ 

^ J_ Ml, 1 

As ^MwTw l+n^ 1+7]'^^ 

- Obtain a (pseudo-)random number pa ; 

- If P3 > P3/P1 map T] — r/. 

- If the new r) falls outside the kinematic limits [77miii, T^max] the event is rejected. 

- Note also that the last mapping can only occur if the original rj was negative, since Ps < Pi only in the region 

Calculate the value of using the inverse of 77 definition: 

sw = [Mw Tw) ■ n + Mw (51) 
The weight corresponding to the sampled value 77 is exactly: 



(52) 



which is the (normalised) inverse of EQuation l48l as requested. 



As it turns out in subsequent generator level studies, this generation procedure gives much better agreement with the differen- 
tial distributions than the usual (width independent) Breit-Wigner; an example obtained for the qq Wbb process is shown 
in Figure|8| The evident consequence is that the unweighting efficiency is substantially improved due to the reduction of the 
event weights in the high s^y region. 
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5.4.2 The Inclusion of Mass Effects in Propagator Sampling 

Studies have shown that the x~" approximation works quite well for t-channel type propagators since the phase space suppres- 
sion factor \/A participates in the denominator, as shown in Eg. 1401 and thus contributes only to the x^" slope. Contrariwise, 
while the a;"" approximation still works reasonably well when describing the s-channel type propagators involving particles 
with high virtuality and/or decay products with low masses, it can be shown that this is not necessarily the case in the LHC 
environment, where the presence of massive decay products can significantly affect the invariant mass distributions. As it can 
be seen in Figure|9|the shape of the propagator dependence can be strongly suppressed by the phase space -\/A (boost) factor 
at low values; thus the sampling function approximation for non-resonant propagators could be approximated with something 
like: 

r , X VA(g,m2,mg) 1 v'A(s,mg,mg) 

/nr(s) = ^ — = (53) 

and similarly 



^ VA(s,mg,mg) ^ \/A(s,mg,mg) 

' S ■ (S - M2)2 + Af2r2 yi-((s-M2)2+M2r2) ^ ' 

for resonant propagators. 

As it turns out the two functions cannot be sampled by the well known unitary algorithms (i.e. the biggest collection of 
recipes |56| yielded no results); already the integral values of the functions yield complicated expressions which cannot be 
easily calculated, let alone inverted analytically. The solution was to code numerical algorithms to calculate the integrals (i.e. 
cumulants) explicitly. 

After the integrals are calculated, their inverse and the subsequent sampling value can again be obtained numerically. 
Namely, resorting to the original definition of the unitarity sampling recipe in Eq. |8| by replacing the normalised gi (x) with: 

9^(^) - ./^""^ . (55) 
/ fix)dx 

which in turn gives: 

X 

f[x) dx = r ■ I f(x) dx^ (56) 



where f(x) is the non-negative function one wants to sample from, [x_,x+] is the range of values of the parameter x we 
want to sample and r a pseudo random number r £ [0, 1]. As already stated (c.f. Eq|8j, in the case when the integral of 
the function f(x) is an analytic function, F{x) = f{x) dx, and has a known inverse F~^(x) one can construct explicit 
unitary prescriptions by: 

X = F-^ (r • [F{x+) - F{x^)] + F{x-)) (57) 

as given for two particular cases in Eq. I42I45I 

In the cases the integral can not be inverted, the prescription of the Eq. I56l can directly be transformed into a zero-finding 
request; thus, since both the integral and the first derivative (i.e. the sampling function and its cumulant) are known, the 
Newton-Rhapson method is chosen as the optimal one for root finding: 



g(x) = { J f{x)dx -r- J f(x)dx } = (58) 



IX + 

ff'(^) = ^ I / .f{x)dx -r- j f{x)dx \ = f{x) 



(59) 



With a sensible choice of starting points the procedure generally takes on the order of ten cycles until finding the root with 
adequate numerical precision. The overall generation speed is still deemed quite reasonable. 
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The integration of the phase-space suppressed resonant propagator of Eg. B^ vields a rather non-trivial expression: 

/r(s) ds = 



\/Hs,ml,ml)ds 



(ma+m5)2 



/J - ((S-M2)2 + M2r2) 



(60) 



\/(s — a)(s — b) ds 
yi- ((s-M2)2 + Af2r2) 



-2ia6r 



-6rM2 (r2 + M2) 

iarcsmh( — ^^j, — 
\/a b 



i arcsinh(- 



—b. a 

7T''b 



+ (i r + M) (a + i (r + i M) M) (6 + i (r + i M) M) n 

+ (r + i M) {b + (-i r - M) M) {ia + {r-iM) M) n 

- (i r + M) (a + 2 (r + i M) M) (6 + i (r + i M) M) n 

- (r + i M) (& + (-1 r - M) M) (i a + (r - i m) m) n 

where the variables a, 6 stand for a — (nia + nib)^ and b — (nia — mb)^ and the functions F[(^, k] and n[<^, k, n] are 
the Legendre's incomplete elliptic integrals of the second and third kind with complex arguments. In order to perform the 
calculations the latter functions had to be coded from scratch since they were not found in any (publicly available) computer 
libraries or code repositories. The prescriptions for calculating them were found in |58 1; the results were checked against the 
values given by Mathematica™. 

In the special case ma = mb the above expression simplifies into: 



M{-iV + M) 
b 


i arcsinh(- 
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'M{ir + M) . 


arcsinh( — 
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/E '■■ 
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M(-ir + M) 
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'M{ir + M) . 
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f r f ^ , f \/Ms,ml,ml)ds 

J ^"(^^^^ ^ J ^.((S-Af2)2 + Af2 



(2ma)2 



(2m„)2 



r2) 



(61) 



\/ {s — a) ds 



•((s- A'/2)2 + M2r2) 

a 

1 

TMy^a + i-iT-M) M 
(i a + (r - j M) M) arctan( 



^a + {-iT~ M) M 



i V'a + (-i r - M) M ^Ja + i{T + iM)M arctan(- 



^a + i {V + iM) M 

The result of integrating the phase-space suppressed non-resonant propagator (Eq.|53j yields a similarly non-trivial result: 



j fNii{s) ds ^ J 



\/Hs,ml,ml) ds 



(62) 



-2 y(a -s)ib- s) F, - (1) , - (1) , 1 - ^, ^, I] 



+ 



V¥^(-ffl + b) {a-s)T\l-u]V[-u,-{\),l-u,l\ 
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where the function F[Qf, 7, x] is the Gauss Hypergeometric function and the Fi [a, f), /3', 7, x, y] is the two-parameter 
(Appell) Hypergeometric function |59|. Both functions can be calculated by using the prescriptions in |59|; it however 
turns out that the calculation of the Fi [a, /3, /3', 7, x, y] to a certain (high) precision is almost two times slower than the 
explicit numerical calculation of the integral to the same precision. Subsequently the numerical evaluation of the Gauss 
Hypergeometric function F[a, /3, 7, x] was retained since it participates in the ma = mb simplification and the calculation of 
the integral was done by using a 50-point Gauss-Legendre quadrature with ^/s weight function; the weights were calculated 
by i60ll . The implementation of the (Appell) Hypergeometric function calculation was used as a cross-check to confirm the 
correct implementation and precision of the numerical method. 
As already mentioned, the above integral again simplifies for ma = mbi 



s 



s 




(63) 



s 




a 



2 





1 5 



and the Gauss Hypergeometric function F[q,/3, 7,x] is in this case calculated by the methods described in |59J with some 
improvements analogous to the ones described e.g. in 1611 . 
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5.5 Application of the Phase Space Generation Algorithms 

The AcerMC 2.0 Monte-Carlo generator uses the multi-channel phase space generation where each channel corresponds to 
an expected phase space topology as derived from the participating Feynman diagrams. In AcerMC 2.0 this information 
was obtained from the modified MadGraph| 8 1 program which also supplied the probability amplitudes for the implemented 
processes. Each channel topology was in turn constructed from the t-type and s-type modules and sampling functions de- 
scribed in this paper together with some additional importance sampling techniques for space angles and rapidity distributions 
described in detail elsewhere 1 1 37 38 1. The unknown slope parameters (denoted v in the text, c.f. Eq. I53t of the invariant 
mass sampling functions for non-resonant propagators were obtained by short training runs of the program on a process by 
process basis. 

As a further step the multi-channel self-optimisation procedure was implemented in order to minimise the variance of the 
event weights further |40|. Eventually, additional smoothing of the phase space was obtained by using a modified VEGAS 
routine to improve the generation efficiency (see next Section). 

The procedure of multi-channel importance sampling used in the event generation can briefly be outlined as follows. The 
analytically integrable function <?($) (c.f. Eq. |9| which aims to approximate the peaking behaviour of the differential cross- 
section dependence on various kinematic quantities is introduced into the differential cross-section equation as a weighted 
sum of several channels (/i(<I'), each adapted to a certain event topology: 

= ^Q,-g,($). (64) 

i 

The values of relative weights ai are determined from multi-channel self-optimisation procedure in order to minimise the 
variance of the weights 1401 . The phase space points are than sampled from the function <?($), first by randomly 

choosing a channel i according to the relative frequencies cxi and then deriving the required four momenta from the chosen 
using unitary" algorithms [371. 
A few representative invariant mass distribution comparisons between the implemented sampling functions and the actual 
differential distributions are shown in Figure|9| 




10 20 30 40 50 60 70 80 90 100 20 40 60 80 100 120 140 160 180 200 20 40 60 80 100 120 140 160 180 200 

m^^ [ GeV/ c- ] [ GeV/ c' ] nij, [ GeV/ ] 



Fig. 9: A few representative invariant mass distribution comparisons between the (normalised) sampling 
functions and the normalised differential cross-section as obtained with AcerMC 2.0 Monte- 
Carlo generator. Left: The invariant mass of the bb pair in the process ud W^g* l^//ibb. 
Center: The invariant mass of the Wbb system (equivalently the hard centre-of-mass energy \/I) 
for the same process. Right: The invariant mass of the pair in the process gg /'y*hb — > 

Uhb. All the distributions were obtained using the prescriptions of this paper without the adaptive 
algorithms also used in the AcerMC 2.0 Monte-Carlo generator. As one can see the approxima- 
tions used seem to work quite well. 



Unitary in this context meaning that there is no event rejection in the algorithm. 
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5.5.1 Modified VEGAS Algorithm 

Using the described multi-channel approach, the total generation (unweighting) efficiency amounts to about 3 — 10% de- 
pending on the complexity of the chosen process. In order to further improve the efficiency, a set of modified VEGAS j?! 
routines was used as a (pseudo-)random number generator for sampling the peaking quantities in each kinematic channel. The 
conversion into a (pseudo-)random number generator consisted of re- writing the calling routines so that instead of passing 
the analysed function to VEGAS for sampling and integration, VEGAS calls produce only (weighted) random numbers in 
the region [0, 1] and the corresponding sampling weight, while the VEGAS grid training is done using a separate set of calls. 

After training all the sampling grids (of dimensions 4-7, depending on the kinematic channel), the generation efficiency 
increased to the order of 6 — 14%. The motivation for this approach was that in unitary algorithms only a very finite set of 
simple sampling functions is available, since the functions have to have simple analytic integrals for which an inverse function 
also exists. Consequently, the non-trivial kinematic distributions can not be adequately described by simple functions at hand 
in the whole sampling domain (e.g. the r distribution, c.f. Figure ITOI and some additional smoothing might be welcome. In 
addition, the random numbers distributions should, due to the applied importance sampling, have a reasonably flat behaviour 
to be approached by an adaptive algorithm such as VEGAS'". 

The further modification of VEGAS, beside adapting it to function as a (pseudo-)random number generator instead of 
the usual integrator, was based on the discussions 1 43 44 1 that in case of event generation, i.e. unweighting of events to the 
weight one, reducing the maximal value of event weights is in principle of higher importance than achieving the minimal 
weight variance. Since the VEGAS algorithm was developed with the latter scope, some modification of the algorithm was 
necessary. As it turned out, the modification was fairly easy to implement: Instead of the usual cumulants: 

< I >cell = ^ Wti, (65) 
coll 

according to the size of which VEGAS decides to split its cells, the values: 

< F >coii = Acoll ■ wt^n" - ^ wti, (66) 

cell 

were collected and used as the splitting criterion. The above value (called loss integral in ]j44jj) is basically a measure of 
the deviation between the maximal weight sampled in the given cell wt^fi" and the average weight in the cell < wtcdi >= 
(X^ccii wti)/Accii (the quantity Acoii denoting the cell width, i.e. the integration range). Re-writing the above expression as: 

<F>ccii = (Acoll -wt-fD-ll-^^^^} (67) 

clearly indicates that the value < F >ccii is actually a measure of the generation inefficiency in the cell, since the term in 
the curly brackets is equivalent to one minus the generation efficiency < wtcoii > /wt™^?. In addition, the inefficiency is 
weighted with the crude/maximal estimation of the function integral over the cell Accii ■ wt™if and cells with the highest 
< F >ccii are split. 

This method is of relevance because the VEGAS cells are actually projections of the whole phase space on the (chosen) 
side axes, i.e. VEGAS cannot isolate a maximal weight in a certain point in phase-space and build a cell around it, which in 
principle would be an ideal solution. An implementation with this scope in view has been made in FOAM 1441 . nevertheless 
we have not found it competitive with respect to the modified VEGAS for the given application. 

The thus modified ac-VEGAS algorithm further increased the unweighting efficiency for almost a factor of two. 

One of the sampling distributions is shown in Fig. 1101 as a gray histogram (marked channel) and the actual (generated) 
differential cross-section dependence is drawn in black. In the first figure, the random variable used for sampling values from 
1/r'' distribution was drawn from a flat probability in the interval [0, 1]; in the second plot the ac-VEGAS algorithm was 
used to give an optimal grid for sampling the random variables needed for parameter generation (the grid is trained for each 
kinematic channel separately, the sum of all channels is shown in the plot). The improvement is evident; one has to stress that 
the use of ac-VEGAS algorithm to generate the values of r directly would be much less efficient since VEGAS gives a grid 
of 50 bins/dimension, which would give a very crude description of the r distribution compared to the one at hand. 

Observing the distributions of the event weights before and after the inclusion of the modified ac-VEGAS algorithm 
(Fig. IIH it is evident that ac-VEGAS quite efficiently clusters the weights at lower values. Note that the principal effect 

'^At this point also a disadvantage of using the adaptive algorithms of the VEGAS type should be stressed, namely that these are burdened with the need 
of training them on usually very large samples of events before committing them to event generation. 
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Fig. 10: Comparison between the sampling distribution for the r = s/ s G [T„iin, 1] variable 
in gg— > ttbb process before and after the application of modified ac-VEGAS STj 
smoothing procedure (light gray histogram). The generated (normalised) differential 
cross-section is also drawn (black histogram, labelled Generatedj. 





Multi-Chan only 




■ Multi-Chane VEGAS 




- Multi-Chaneac-VEGAS 







Wtx5 • lO' 



Fig. 11: The distribution of event weights using only the Multi-Channel approach (dotted his- 
togram) and after application o/ VEGAS (dashed histogram) and ac-VEGAS (full his- 
togram) algorithms in the gg (Z*^ ~*)llhh process. 



of original VEGAS is indeed to cluster event weights in a narrow region, nevertheless a tail towards the high-weight region 
remains. On the other hand, the ac-VEGAS efficiently reduces the tail in the high weight region; only a few of the event 
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weights still retain their large values, thus reducing the generation efficiency. Given the difference in distributions, the ob- 
served increase of the generation efficiency seems relatively modest. To better understand this result one should consider that 
the formula for the MC generation efficiency is given by: 

< wt > 

e = — , (OS) 

Wtmax 

where < wt > is the average weight of the sample and equals the total event cross-section, while wtmax represents the 
maximum event weight in the applied generation procedure and is determined through a pre-sampling run with a high statistic. 
Since the average weight < wt > equals the total cross-section of the process, it remains (necessarily) unchanged after the 
application of the VEGAS refining; consequently the change of efficiency results in the reduction of the maximum weight 
wtmax by approximately a factor two, which is from technical point of view quite an achievement. 

A further step to profit from the clustering of weights induced by ac-VEGAS is to adopt a re-definition of the MC 
generation efficiency as proposed by 1421 1441 . In this approach, the alternative definition of wtmax is: For a given precision 
level a « 1, the wtmax is determined from the total weight distribution in such a way that the contribution of the events 
exceeding this value to the total weight sum (i.e. cross-section integral) equals a. Such a quantity is referred to as wtmax ™d 
the efficiency expression becomes: 

(69) 



The argument presented in f42"'441 seems to be quite reasonable since the true event weight is in any case only estimated 
from a finite sample of events and the new definition simply takes into account a certain level of accuracy in the maximum 
weight determination. In addition, certain very weak singularities that might exist in the simulated process and might occa- 
sionally result in a very high event weight are automatically taken into account. The use of new wtmax consequently results 
in a generation efficiency of about e > 20% for all the implemented processes, which is a significant improvement in terms 
of time needed for MC generation. 



5.5.2 Colour Flow Information 

Before the generated events are passed to PYTHIA/HERWIG to complete the event generation, additional information on 
the colour flow/connection of the event has to be defined. Below we discuss the implemented method of the colour flow 
determination on the example of two processes, gg — > tibb and qq — > tibb. 

For the process gg tibb six colour flow configurations are possible, as shown in FigureQ2| With 36 Feynman diagrams 
contributing to the process and at least half of them participating in two or more colour flow configurations, calculations by 
hand would prove to be very tedious. Consequently, a slightly modified colour matrix summation procedure from MAD- 
GRAPH 1 8 1 was used to determine the colour flow combinations of the diagrams and the corresponding colour factors. The 
thus derived squared matrix elements for separate colour flow combinations |A^flowl^ were used as sampling weights on an 
event-by-event basis to decide on a colour flow configuration of the event before passing it on to PYTHIA/HERWIG for 
showering and fragmentation. The procedure was verified to give identical results regarding the colour flow combinations 
and corresponding colour factors when applied to the processes published in |45 1. As one can see this approach neglects the 
interference terms between the distinct colour-ordered amplitudes and is indeed exact only in the Nc — > oo limit'^^ I46I I47I . 
Since there is no a priori rule of how to split the interference terms between the colour ordered amplitudes this approach is 
generally deemed to be the best one can do; recent developments in this field 1481 however suggest additional improvements 
to the method that indeed might be incorporated into later versions of AcerMC. 

The colour flow configuration in the qq tibb channel is much simpler since only two colour flow topologies exist (Fig. 
I13t ; the choice between the two has been solved in a manner identical to the one for the gg — > tibb process, as described 
above. 

Some specifications of the implemented matrix-element-based processes: number of Feynman diagrams, channels used 
in the phase-space generation and colour flow configurations are collected in Table ll2l The control processes ID=91 -94 are 
omitted from the table due to their simplicity. 



5.6 The ckqed and cxs calculations 

Native functions of running qqed(Q^) and as(Q^) have been implemented inside AcerMC with the main objective of 
providing a means to keep the (total) cross-sections of the processes unchanged when interfacing with the two supervising 

The matrix elements used in the cross-section calculation and event generation are of course complete and do not employ any approximation. 
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Fig. 13: A diagrammatic representation of the two colour flow conflgurations in the process 
qq ttbb. 



generators, since the implementations of the two functions in PYTHIA and HERWIG differ to some extent. Especially the 
as{Q^) is subject to experimental and theoretical uncertainties, however obtaining a different cross-sections for the same 
AcerMC process due to different interface, could be regarded (at least to some extent) as an inconsistency 

• oqed is implemented in AcerMC using the formulae given in 1 62 1 and is in complete accordance with the implemen- 
tations in PYTHIA and HERWIG apart from the updated hadronic component published recently by Burkhardt et. al. 
|63). As one can see in Figure [l4l the latter minimally lowers the qqed values. 

• a.s has one and three loop implementations in AcerMC following the calculations of W. J. Marciano (SM and using 
A^jg' transformations for flavour threshold matches. The three loop version gives good agreement with the HERWIG 
implementation (both functions have been set to the same A<%=^) value) as one can see in Figure lT31 The PYTHIA 
two loop implementation deviates somewhat from the latter two; the kinks observed in the plot are due to approximate 
^aJs' transformations at flavour thresholds, which are exact to one loop only. 

Although the AcerMC and PYTHIA one loop implementations are identical in form the resulting values differ by a small 
amount because the default PYTHIA implementation reads the A^^*^^^^' value from PDFLIB instead of the Aj^^^^' one used 
by AcerMC and HERWIG; the difference thus occurs due to A^'^^' propagation at flavour thresholds. 

'"'The values will still differ by a small amount in processes containing W bosons (processes 3,4) due to different values of the CKM matiix in the two 
supervising generators. 
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Table 12: Some details on matrix-element-hased process implementation in AcerMC library. 

In case ofqq initial state the number ofFeynman diagrams corresponds to one flavour 
combination. The f = e,iJ,, r, b. 
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Fig. 14: Comparison between the aQED(Q'^) fLeftj and as{Q^) (Rightj implementations in 
AcerMC, PYTHIA and HERWIG. For as(Q^) calculations with different loop orders 
(L) are given where applicable. 
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6 Structure of the package 

The AcerMC package consists of a library of the matrix-element-based generators for selected processes, interfaces to the 
PYTHIA 6.2 and HERWIG 6.5 generators, sets of data files and two main programs: demo_hw.f and demo_py.f. Provided 
makefiles allow to build the executables with either of these generators as the supervising generator, demo_hw.exe and 
demo_py.exe. 

6.1 Main event loop and interface to pythia/herwig, tauola and photos 

The main event loop is coded in the demo_hw.f or demo_py.f files, where the opening/closing of the input/output files, 
reading of the data-cards and event-loop execution is performed. Main event loop consists only of calls to the acermc.py or 
acermC-hw subroutines, with parameter MODE = -1,0,1 respectively set for initialisation, generation and finalisation of 
the event loop. The call to acermc_XX activates respective procedures of the supervising generator, which in turn activates 
the acevtgen procedure steering the native AcerMC generation of the matrix element event. Fig.^|illustrates this calling 
sequence in some details. 

As one can deduce from the diagram in Fig. 1151 certain functions called by AcerMC, as e.g. pseudo-random number 
generator acr are re-routed through the interfaces to the linked supervising generator, depending on the choice at compilation 
time (e.g. acr function giving (pseudo-)random numbers is linked to either pyr or hwrgen as shown in the plot), providing 
the internal consistency of the package. The generated event is rewritten to the format required by the supervising generator 
by means of the acdump_XX routines. 

The pythia_ac.f and herwig_ac.f files contain sets of re-routing/interface functions, specialised for the respective supervis- 
ing generator. The main library of AcerMC is well screened from dependencies on the supervising generator, all dependencies 
are hidden in herwig_ac.f and pythia_ac.f respectively. 

The PYTHIA and HERWIG libraries remain essentially untouched'^, without introducing any dependencies on the Ac- 
erMC code. The input cards are common for both interfaces. 

The calling sequence is further enhanced by the optional calls to the external TAUOLA 1 5 1 library, which handles the 
r-lepton decays using the spin information of the hard process, and the PHOTOS |4| library which adds the final state QED 
radiation to final state leptons and hadrons. The two interfaces are controlled by two additional input cards. The details of the 
interface are given in Section FOl 



'^For specification of exceptions see Section 7.4 and 7.5. 
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Fig. 15: The calling sequence of the main event generation routine acevent_xx. The routine is called either through demo_py 
acermc_py sequence when interfacing the PYTHIA 6.2 generator or demo_hw acermc_hw sequence when the HERWIG 6.5 
is linked. The structure of the interface subroutines and relations with the corresponding ones from supervising generators 
and/or external libraries is also evident. 



6.2 Structure of the AcerMC matrix-element and phase-space code 

The AcerMC core code performs the generation of a matrix-element-based event. Fig.Q6|illustrates the caHing sequence for 
generating an unweighted event. The steering subroutine is called acevent_gen which subsequently constructs the weight by 
calling the Madgraph/HELAS subroutines for the matrix element evaluati To stress again, this subroutine calls only a sequence 
of the native AcerMC subroutines, any call to the supervising generator goes via the respective interface function/subroutine. 
A more detailed representation of calling sequence is shown in the Figure [l6l 

Code for the phase space generation is since the version 2.0 greatly simplified compared to the earlier versions. Code for 
matrix element calculations is grouped together for all processes in subdirectory matel. Code with different utility subroutines, 
e.g. kinematic transformations used by all subprocesses, is in the subdirectory connmon. Subdirectory interface contains 
code with interfaces to supervising generators, finally subdirectory include contains all include files. The overall view on the 
structure of the AcerMC directories is shown in Fig. 1171 

The core code builds one library ilbacermc.a. 
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ac-VEGAS 
grids 



Fig. 16: The event generation sequence controlled by acevent_gen subroutine. Phase space generation is sequenced by 
calling the ackinch_gen2 routine to obtain the incoming gluon and the outgoing four-momenta of the participating particles. 
The latter routine handles the possible momenta permutations and calls the explicit four-momenta generation (and PS weight 
calculation) in the sequence prescribed by the event topology. These { channel-specific) routines are constructed from common 
building blocks listed in the next two columns. The acevent.gen routine also initialises MADGRAPH/HELAS package and 
retrieves the matrix element values. All the generated four-momenta, as well as the event weight are finally passed back to the 
supervising generator via the acevent_dump call. 



6.3 Data files for the phase-space optimisation 

The AcerMC matrix-element-based generators are very highly optimised, using multi-channel optimisation and additional 
improvement with the ac-VEGAS grid. The generation modules require three kinds of the input data to perform the generation 
of unweighted events: 

• A file descrbing the construction of the topologies relevant for the chosen process, the implemented sets are stored in 
the directory acermc.dat/tops. 

• A file containing the list of the values of relative channel weights obtained by the multi-channel optimisation, defaults 
being stored in acermc.dat/m-cwt. 

• A file containing the pre-trained ac-VEGAS grid, the pre-trained (default) ones located in acermc_dat/grids. 
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Fig. 17: The structure of the AcerMC directories. 



• A file containing the maximum weight wtmax, a-cutoff maximum weight wt^^x and the 100 events with the highest 
weights, the default ones being provided in acermc_dat/maxims. 

In case of changing the default running conditions, like parton density functions or centre-of-mass energy, the user should 
repeat the process of preparation of the listed data files containing the inputs for the phase-space generator modules in order 
to preserve the initial event generation efficiency. 

The reading sequence of data files inside AcerMC is shown in Figure llSl 
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Fig. 18: The reading sequence of the input files and the performing subroutines in the AcerMC code. 



Pre-trained data sets are obtained using = 14 TeV, PYTHIAdefault as (Q^) and aQED((5^) and CTEQ5L (parametrised) 
parton density function set and are provided for each implemented process'^. For these, the relative channel weights are stored 



"These can also be used for a series of other settings, see Section l7!4l for details 
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in the INCLUDE files in acermcsrc/include/chanwt-XX.inc where XX denotes the process id (c.f. TablelTIl: the default/pre- 
trained ac-VEGAS grids are listed in the directory acermc_dat/grids/vscalA_XxYYY.veg, where A denotes the scale choice 
of the process XX and YYY denotes the cutoff value of the m^o for the AcerMC processes XX = 05 ^ 08. The files con- 
taining the maximal weights wtmax and wt° as well as the 100 events with the highest weights are stored in the directory 
acermc.dat/grids/vtmaxA_XxYYY.dat, following the same labelling convention. Both the trained ac-VEGAS grids and the 
weight files were obtained from test runs with at least 2 • lO" weighted events being generated. 

The number of required input files might at first look seem large, considering that many event generators do not require 
any input files for operation; the difference is not in so much in the complexity of the phase space generation as in the fact 
that many event generators require a warming run instead, i.e. before the generation of unweighted events is performed a 
certain number of weighted events (typically of the order of IC*) is generated in order to obtain the relative multi-channel 
weights (in case multi-channel phase space generation is used) and/or the optimised VEGAS grid and/or an estimate of the 
maximal weight. Such an approach can have an advantage when event generation is very fast and the phase space regions 
with the highest weights are well known (as done for the 2^2 processes in PYTHIA); on the other hand, when the phase 
space topology of the process is more complex and the event generation is comparatively slow, generating a relatively small 
number of e.g. 10* weighted events every time a generator is started can become CPU wasteful and/or inaccurate in terms of 
maximum weight estimation. 

Reasonably accurate estimation of the latter is namely crucial for correct event unweighting; event generators using 
warming-up method for maximal weight search often find still higher weights during the production run and reset the maximal 
weight accordingly. In this case however, statistically correct approach would be to reject all events generated beforehand and 
start the event generation anew, which is almost never implemented due to the CPU consumption and the possibility of hitting 
a weak singularity (the same argument leads to the definition of the wtJJiax. c.f. Section l5.5.1> . With a small pre-sampled set 
the generator can however badly under-estimate the maximum weight and a large number of events can be accepted with a 
too-high probability. The only hope of obtaining correct results is in such cases that the weight plateau will be hit sufficiently 
early in the event generation process. Consequently, such approach can be very dangerous when generating small numbers of 
events". 

In contrast to the warming-up approach for the AcerMC we decided that using separate training runs with large numbers 
of weighted events to obtain the optimised grids and maximum weight estimates are preferable. In case a user wants to produce 
data sets for non-defaults setting, this can easily be done by configuring the switches in the acermc.card (see Section 6.2). 



Small being a somewhat relative quantifier, since the size of an representative sample should depend on the phase space dimension, i.e. the number of 
particles in the final state; with e.g. 4 particles in the final state, 10^ events can still be considered a relatively small statistics. 
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7 How to use the package 

There are two main steering input files: run.card and acermc.card which share a common format for both executables. The 
run.card (see App. IB.H provides switches for modifying: generated process, number of events, parton density functions, 
predefined option for hadronisation/fragmentation in the supervising generator, random number, etc.. The acermc.card (see 
App. \B.2i provides switches for modifying more specialised settings for the AcerMC library itself. Once the user decides 
on a setup for the generated process, only run.card is very likely to be modified for the job submission. In case the user 
decides to use the external TAUOLA or PHOTOS libraries by selecting the appropriate switches in the run.card there are 
two additional files tauola.card(see App. IB.3> and photos.card (see App. IB.4t which steer the performance of these two 
libraries. All input files are read by AcerMC executables through the CERNLIB FFREAD routines, some commands given 
in the input files (e.g. LIST entry, see Appendix B) are internal FFREAD commands which should be disregarded by the user. 

The same executables can also be used for running standard PYTHIA 6.2 and HERWIG 6.5 processes. The example how 
to require such process is provided as well, in demo.hw.f and demo_py.f respectively. If the user requires that the AcerMC 
library is not used, the ttH production will be generated with demo_py.exe and HERWIG 6.5 implementation of the Zbb 
production will be generated with demo_hw.exe. In this case only the run.card file will be read, so in case the user requires 
different processes and/or settings of the supervising generators the user has to implement her/his steering there or create 
another xxx.card file, together with the corresponding code added to the demo_XX.f. 



7.1 Steering switches of the overall run 

The overall run is controlled by the switches read from the run.card file (see also Add. IB.lV Some of these general switches 
are also passed to the AcerMC library. 

. CMS : Sets the centre-of-mass energy in GeV. 

. ACER : Specifies if the internal AcerMC process will be used 
ACER=0 - use process from PYTHIA/HERWIG 
ACER=1 - use internal AcerMC process 

ACER=2 - use internal AcerMC process & dump events into a record 
ACER=3 - use internal AcerMC process & read events from a record 

. PROCESS: Sets process id 



• HAD : Sets predefined option for QCD ISR/FSR and hadronisation 
HAD=0-only hard process 

HAD=1 - only ISR (works for PYTHIA interface only) 

HAD=2 - only ISR and FSR 

HAD=3 - full treatment 

HAD=4 - only FSR 

HAD=5 - only FSR and hadronisation 

. PDFLIB804NGROUP 

Sets the value of the PDFLIB804NGROUP parton density function choice. 

. PDFLIB804NSET 

Sets the value of the PDFLIB804NSET parton density function choice. 
. RSEED : Choose the random seed for (pseudo-)random generator initialisation 



. TAUOLA: Specifies if the TAUOLA library will be used for tau decays 

TAUOLA=0 - use internal PYTHIA/HERWIG mechanisms for r-lepton decays 
TAU0LA=1 - use the TAUOLA library for rJepton decays. 



. PHOTOS : Specifies if the PHOTOS routines will be used for final state photon radiation 
PHOTOS=0 - use internal PYTHIA/HERWIG mechanisms for FSR photon radiation 
PH0T0S=1 - use the PHOTOS mechanisms for FSR photon radiation. 



• N EVENT : Required number of generated events 
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7.2 Steering switches of the AcerMC processes 

The AcerMC processes are controlled by values set in a simple arrays specified in acermcsrc/include/AcerMCinc: 
C CROSS-TALK PARAMETERS 
DOUBLE PRECISION ACSET 
INTEGER lACPROC 

COMMON/ACPAR1/ACSET(200)JACPROC(200) 

C PARTICLE PROPERTIES 

DOUBLE PRECISION ACCHG,ACMAS,ACCKM 

COMMON/ACPAR2/ACCHG(50,4),ACMAS(50,4),ACCKM(4,4) 

C ROUTINE I/O 

INTEGER LACSTD,LACIO 

COM MON/AC PAR3/LACSTD, LACIO 

The lACPROC array activates the process IPROC=PROCESS (read from run.card file) by setting IACPR0C(IPR0C)=1 . 

The list of currently implemented processes in AcerMC can be found in Table llll When running in the generation mode 
with ACER=0 full list of processes implemented in either PYTHIA 6.2 or HERWIG 6.5 can be activated, however the 
mechanism for passing information about process id to either of these generators has to be coded by user individually in 
demo_xx.f. 

The main control switches reside in the array ACSET. The COMMON block ACPAR2 contains the particle charges, 
masses and decay widths as well as the CKM matrix using the PYTHIA convention. The values are filled by the interface 
routines to be equal to the PYTHIA/HERWIG internal values in order to preserve consistency within the generation stream. 
In case the user wants to change some of the particle properties this should be done through the native PYTHIA/HERWIG 
switches; AcerMC will copy them and use the new values. 

The COMMON block ACPAR3 contains the two logical I/O unit numbers used by AcerMC. The LACSTD value deter- 
mines the output unit of the AcerMC messages and the LACIO unit is used for readi ng/writing the AcerMC data files. 

The main control switches which reside in the array ACSET (see also App.IR21: 

. ACSET(1) : Sets the centre-of-mass energy in GeV. 

. ACSET(2): Scale of the hard process 

Choose the sca le for the active AcerMC process. The implemented values differ for various processes, the currently implemented settings are speci- 
fied in Section lT^ 

. ACSET(3): Fermion code 

The flavour of the final state fermions produced in . /'j* // decays of AcerMC processes 3 — » 8. The PYTHIA/PDG naming conven- 
tion is used: 

ACSET(3) = 1 - W ^qq- Z°/7* ^ qq 
ACSET(3)=4-H' ^ ei/,,pi.,,;Z°/7* ^ e+e-,M+M- 
ACSET(3) = 1 1 - W ^ ei.,; Z«/7* e+e- 

ACSET(3)=12- 

ACSET(3)=13-iv-»p''m;^V7* ^ p+p" 

ACSET(3) = 1 5 - IV ^ Ti.^; Z'V7' ^ r + r- 

ACSET(3)=5-zV7- ^ bb 

At present the ACSET(3)=5 is implemented only for processes 7—^-8. 

. ACSET(4): z"/7* propagator 

Use full propagator instead of the pure propagator in matrix element calculation for the AcerMC processes 5 — > 8 and 91. The switch 

is provided since in some of the analyses the 7* contribution is of relevance in the selected mass windows; for the analyses selecting the mass window 
around the Z^ peak this contribufion can safely be neglected. 
ACSET(4)=0 - only Z" propagator. 
ACSET(4) = 1 -fufl ZV7* 

propagator. 

. ACSET(5): 

^zO/-y* mass cut 

Cutoff value on the invariant mass "T'^Oy-y* in GeV when ACSET(4)=1 . Note that the provided data files exist only for values of ACS ET(5)=2,5, 10, 15,30,60, 
120,270,300 and 500 GeV which should satisfy most user requirements for the analyses foreseen at LHC. In case a different value is set the user has also 
to provide the user data files for the mn. 

• ACSET(6) : Sets the value of the PDFLIB804 NGROUP parton density function choice. 

. ACSET(7): Sets the value of the PDFLIB804NSET parton density funcdon choice. 
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• ACSET(8): The implementation of a, (Q^) 

Selects the implementation of as (Q^ ) to be used in the matrix element calculation: 
ACSET(8)=0 - Use the (Q^) as provided by the supervising generator 

ACSET(8)=1 - Use the Ois (Q^) (one loop) provided by the AcerMC; this option gives (Q^) values equal to the default PYTHIA implemen- 
tation. 

ACSET(8)=2 - Use tiie as (Q^) (thi«e loop) provided by tiie AcerMC. 

. ACSET(9):a,(M|)value 

Sets the CKs (M^) value to be used in the ex., (Q^) calculations in case the AcerMC native implementation (ACSET(8)=1 ) is used. 
ACSET{9)=-1 - The value is taken from the PDFLIB804 for the selected parton density function set. 

ACSET(9)>0 - The provided value is taken. 

• ACSET{1 0): The implementation of qqed (Q^) 

Selects the implementation of aqED (Q^ ) to be used in the matrix element calculation: 
ACSET{1 0)=0 - Use the oqedCQ^) as provided by the supervising generator. 
ACSET(10)=1 -Use die aQED(Q^) implemented in die AcerMC. 

. ACSET(11) : c«qed(0) value 

Specifies the value of oqed (0) for AcerMC oqed (Q^) calculation. 
ACSET(11)=-1 -Theo!QED(0) valueis setto o;qed(0) = 0.0072993. 
ACSET(1 1 )>0 - The provided value is taken. 

. ACSET{12) : Decay mode of the produced tt pair 

Sets the decay mode of the W boson pair from the if final state in the AcerMC processes 1,2,4,7,8,9,92 and 93. For ACSET(1 2)>0 the combinatoric 
value of the a X Si? is recalculated and printed in the output. This switch was implemented since the supervising generators (PYTHIA/HERWIG) 
do not allow for forcing specific decays of the top quark pairs generated by external processes. This switch imposes a modification of the decay tables of 
the supervising generators on an event by event basis. 

ACSET{1 2)=0 - both W bosons decay according to PYTHIA/HERWIG switches. 

ACSET{12) = 1 - Wi eur_ and W2 ^ qq. 

ACSET{1 2)=2 - Wi and W2 ^ qq. 

ACSET{1 2)=3 - Wi Tv^ and qq. 

ACSET{12)=4 - eur_,^iy^ and W2 ^ qq. 

ACSET{12)=5 - eiy^,fiiJf, andVl^2 ^ qq- 

The settins ACSET{12)=5 works for PR0CESS=4 only and implies leptonic decay of the 1^ -boson with the same charge as the one of the 
primary W boson produced in the hard process. Folowing configurations are posible: 
qq W+ti^ {W+ -*)L+vl (W+ -*)l+vi b (W~ -^)q'q' b, 
or: 

qq ^W^tl^ {W- ^)L-iJL (W+ ^)q'q' b (M'7^)rp, b, 

where is the lepton from the primary W decay (controlled by ACSET(3) switch) and l'^ is either an or as for ACSET(1 2)=4. Since 
the charge of the semi-leptonic decaying W is correlated with the charge of the primary W boson, the (t X Si? is consequentiy a factor two smaller than 
the oneforACSET{12)=4. 

• ACSET(1 3): Decay mode of the produced WW pair, works for PROC=ll-14 only 
ACSET{1 3)=0 - both W bosons decay according to PYTHIA/HERWIG switches. 
ACSET{1 3)=1 - Wi — ei^e and W2 qq. 

ACSET{1 3)=2 - Wi iJ.Vi_, and W2 -> qq. 
ACSET{1 3)=3 - Wi Ti/^ and W2 qq. 
ACSET{13)=4 Wi ^ ei/e./ii^^ andW2 ^ qq- 

ACSET{1 3)=1 1 -bofli w^eue. 
ACSET(13)=13-boaiiy ^ /m.^. 
ACSET(13)=15-boaiiy 

ACSET(1 3)=20 -bofll W ^ qq. 

• ACSET(50): AcerMC ttaining mode 

The switch controls the mode in which AcerMC is run: 
ACSET(50)=0 - production run, generate unweighted events. 

ACSET(50)=1 - perform multi-channel optinusation and output the user file with channel weights. 
ACSET(50)=2 - perform aC-VEGAS grid ttaining and output the user file with u^ed aC-VEGAS grid. 
ACSET(50)=3 - perform aC-VEGAS grid ttaining as in ACSET(50)=2 but do this by updating a provided grid. 

• ACSET(51 ): Required number of generated events NEVENT 

In case the switch ACSET{50) is set to the non-zero value (i.e. in one of the training modes) the ACSET(51 ) entry is used and defines the number 
of (weighted) events that will be generated; this information is necessary for the learning algorithms to decide on steps in the learning sequence. 
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• ACSET(52): User data files 
Use the data files provided by user: 

ACSET{52)=0 - no, use native (default) AcerMC data files. 

ACSET(52)=1 - use the user's multi-channel optimisation and VEGAS grid files. 

ACSET{52)=2 - use the default multi-channel optimisation and user's VEGAS grid files. 

ACSET{52)=3 - use the default multi-channel optimisation and VEGAS grid files; read the user maximal weight file. 

. ACSET(53) : Maximum weight search 

Mode for the maximum weight search needed for unweighting procedure: 
ACSET(53)=0 - no, use the provided files containing maximal weights. 

ACSET{53) = 1 - use the provided files for max. weights and re-calculate the max. weights using the stored 100 events with the highest weight. 
ACSET(53)=2 - perform the search and give the new Wtm3X_XX_nGW.d3.t file; the switch is equivalent to generation of weighted events. 

. ACSET(54) : Maximum weight choice 

Use the a-cutoff maximal weight wt^^^, or the overall maximal weight w^tmax found in training (see Section |5.5. Il for the explanation on these two 
options). 

ACSET(54)=0 - use the wt" weight. 
ACSET(54)=1 - use the wt^ax weight. 

• ACSET{56): Naive QCD conection for width calculations 

Use the naive QCD multiplicative corrections in the resonant width calculations of the weak boson and top quark decay widths. It is recommended to use 
them except for some specific cases or studies. 
ACSET(56)=1 - use the naive QCD conections 
ACSET(56)=0 - don't use naive QCD conections 

7.3 Steering tauola and photos 

It is highly recommended that in case TAUOLA is used for r decays that PHOTOS is also activated and the setting PMODE 
is set at least to 2 in order to keep the r-leptons radiating (see below), since it gives the necessary contribution to the expected 
T branching ratios. 

It should also be stressed that the actions of the interfaced TAUOLA and PHOTOS are meaningful only when the full 
hadronisation procedure of the event is selected (by setting HAD 3 in run. card). 

The switches and settings specific to the TAUOLA library are set in the tauola.card: 

. POLAR : Polarisation switch for tau decays 
POLAR=0 - switch polarisadon off 
P0LAR=1 - switch polarisadon on 

. RADCOR : Order(alpha) radiative con'ections for tau decays 
RADCOR=0 - switch coiTcctions off 
RADC0R=1 - switch conections on 

. PHOX : Radiative cutoff used in tau decays 

PHOX=0.01 - default value by TAUOLA authors 

. DMODE : Tau and tau pair decay mode 
DMODE=0 - all decay modes allowed 

DM0DE=1 ' (LEPTON-LEPTON): only leptonic decay modes 
DM0DE=2 - (HADRON-HADRON): only hadronic decay mode 

DM0DE=3 - (LEPTON-HADRON): one tau decays into leptons and the other one into hadrons 
DM0DE=4 - (r — i- -Kv) : taus are restricted to decay to a pion and neutrino 

. JAK1/JAK2: Decay modes of taus according to charge, the list is taken from the TAUOLA output. The listing gives only r modes, the are 
charge conjugate, neutrinos are omitted. JAK1 : decay mode of r~^, JAK2 : decay mode of r~ : 



JAK1/2 = 


1 - 


T~ 


e— 


JAK1/2 = 


2 


T~ 


M- 


JAK1/2 = 


3 


T~ 


^ 7T 


JAK1/2 = 


4- 


T~ 




JAK1/2 = 


5- 


T~ 


— > (two subch) 


JAK1/2 = 


6 


T~ 




JAK1/2 = 


7- 


T~ 


— > K*^ (two subch) 


JAK1/2 = 


8 


T~ 




JAK1/2 = 


9 


T~ 




JAK1/2 = 


10 


' - r 
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JAK1/2 




11 - 


T~ 


— * 


37T- ,2n+ 


JAK1/2 




12 


r~ 




37r~, 27r+, tt" 


JAK1/2 




13 


T~ 




27r~, 7r+, 37r° 


JAK1/2 




14 


r~ 




K~,7r~,K+ 


JAK1/2 




15- 


T~ 






JAK1/2 




16- 


r~ 






JAK1/2 




IT- 


r~ 






JAK1/2 




IS 


T~ 






JAK1/2 




19- 


T~ 






JAK1/2 




20- 


T~ 






JAK1/2 




21 - 


T~ 






JAK1/2 




22- 


T~ 




K-K° 



The switches and settings specific to the PHOTOS routines are set in the photOS.Card: 

. PMODE : Radiation mode of photos: 

PM0DE=1 - enable radiation of photons for leptons and hadrons 
PM0DE=2 - enable radiation of photons for r-leptons only 
PM0DE=3 - enable radiation of photons for leptons only 

. XPHCUT: Infrared cutoff for photon radiation: 

XPHCUT=0.01 - default value by PHOTOS authors 

• ALPHA: Alpha(QED) value used in PHOTOS: 
ALPHA j - leave default (0.00729735039) 

• INTERF : Photon interference weight switch 
INTERF = 1 - interference is switched on 
INTERF = - interference is switched off 

. ISEC : Double bremsstrahlung switch: 

ISEC=1 - double bremsstrahlung is switched on 
ISEC=0 - double bremsstrahlung is switched off 

. IFTOP: Switch for pp(gg) ft process radiation: 
IFT0P=1 - the procedure is is switched on 
IFTOP=0 - the procedure is is switched off 



Detailed information about TAUOLA and PHOTOS implemetations in ttie AcerMC setup can be found in Sections [7. 101 and 

ItTTI 
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7.4 How to prepare data-files for the non-default setup 

The following actions are possible, to recover better efficiency of the generator modules with the non-default settings: 

• The user wants to generate events using different parton density function sets and/or different coupling values (e.g. 
AcerMC third order Us {Q^) instead of the first order one): 

It should suffice to set the the switch ACSET(53)=1 , which signals AcerMC to re-calculate the wtmax and wt™ ax 
using the 100 events stored in the file 

acermc_dat/grids/vtmaxA_XxYYY.dat. The coupling and parton density functions values should not change signifi- 
cantly the process topology but affect foremost the overall scale of the event weights; thus, the stored hundred events 
should still remain the ones with the highest weights and the re-calculated approximate estimates of the highest weight 
should be accurate enough. 

In case the user is not confident in the obtained result, the new maximal weight estimation can be initiated by setting 
the switch ACSET(53)=2, which will result in generation of weighted events. The number of generated events is deter- 
mined by the usual NEVENT in run. card. At the end of the run AcerMC will produce a file called wtnnax_xx_new.dat, 
with XX specifjdng the process number. The user should then start the generation of unweighted events with the setting 
ACSET(52)=3 and linking(renaming) the new file to wtmax_xx_usr.dat, with xx denoting the process number (e.g. 
wtmax_01_usr.dat). 

• The user wants to generate events using different values of particle/boson masses or other significant changes of the 
parameters apart from the centre-of-mass energy and/or m zo cutoff value for processes 5-8: 

In this case the user should re-train the VEGAS grid since the process topology is assumed to undergo minor changes. 
This is done by setting the switch AGSET(50)=2 or ACSET(50)=3; in the first case AcerMC starts with an untrained 
grid and in the second one it starts modifying the existing grid provided for the process at the selected hard process 
scale. In general the second option should be preferable since the topology should still be close to the pre-trained one. 
AcerMC again produces weighted events and at the end of the run outputs a file grid_xx_new.veg. The number of 
generated events is determined by the usual NEVENT in run. card. As in the previous case, the user should re-name the 
file to grid_xx_usr.veg, re-set the switch to ACSET(50)=0 and repeat the maximal weight search procedure described 
above, by setting the switch ACSET(53)=2 etc.. When the maximum weight search is completed the user switch AC- 
SET(52)=2, which will cause AcerMC to read the wtmax_xx_usr.dat as well as grid_xx_usr.veg files and produce 
unweighted events with the new setup. 

• The user wants to generate events at a different ruzo /.y* cutoff value and/or different centre-of-mass energy y^: 

When the user changes at least one of these two parameters the event topology is significantly changed as well as the 
contributions from different kinematic channels. The user should thus start with a new multi-channel optimisation by 
setting the mode switch AGSET(50)=1 and start an AcerMC run. The number of generated events is determined by the 
usual NEVENT in run. card. At the end of the run AcerMC will produce a file chanwt_XX_new.dat which should be 
renamed/Unked to chanwt_xx_usr.dat. The user should then set the switch AGSET(52)=1 and first put AGSET(50)=2 
and and repeat the VEGAS grid training as described above and consequently ACSET(50)=0 and AGSET(53)=2 to 
perform the maximum weight search. After obtaining all three user files the AGSET(53) should again be put back to 
ACSET(53)=0 and a normal run should be started; the switch ACSET(52)=1 will in this case force AcerMC to read 
all three user files and produce unweighted events. 

At the first look procedure for listed action scenaria might seem a bit complex but should after a few trials and errors become 
a straightforward routine; it is expected that the vast majority of users would have to deal with at most the first scenario. 
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7.5 Details on the interface to PYTHIA 6.2 

The AcerMC interface to Pythia 6.2 is implemented close to the new standard specified at the Les Houches workshop 
2001 |3|. The full description of the standard can be found in the PYTHIA 6.2 manual (| 1 1). In addition to the UPINIT 
and UPEVNT routines the file acermcsrc/interface/pythia^cf provides links between a list of AcerMC routines and the 
corresponding PYTHIA ones, as e.g. the (pseudo-)random number generator, as and ciqed calculations as well as a series 
of routines that re- write the AcerMC event output to the required PYTHIA format. Using this strategy, the native AcerMC 
code is completely de-coupled from the linked hadronisation library (at the moment PYTHIA/HERWIG) and new interfaces 
can thus easily be added. The special AcerMC requirement is the call to the AGFINAL subroutine at the end of the run which 
signals the AcerMC to close the various I/O files and produce the final output. An example of the implementation of the 
PYTHIA/AcerMC interface can be found in the provided demo_py.f. The PYTHIA code is unmodified apart from making a 
small modification in PYINIT routine: 
CALL UPINIT(I) 
..parameter initalisation.. 
CALL UPINIT(2) 
..process initialisation.. 

since the user-supplied processes in this new interface are not allowed to (re-)estimate maximal weights (as e.g. the native 
PYTHIA processes do). In the original code the call to UPINIT is set before the PYTHIA parameters and functions (e.g. 
PYALPS for Qs(Q^) calculation) are initialised with the user settings . 

An additional modification of the PYTHIA code was needed in order to put the top decay products (handled by PYTHIA) 
back into the history part of the PYTHIA event record, which was deemed useful for e.g. top reconstruction studies. The 
modification is however assumed to be only temporary, i.e. until the feature is added in the PYTHIA code itself. 

The user can thus add the most recent PYTHIA library without other neccessary modifications but for the two lines of 
code in PYINIT routine as described above (the dummy routines UPINIT, UPEVNT, STRUCTIVI, STRUCTP and PDFSET 
however have to be removed from the code for an external process to work and to activate the PDFLIB interface). 

By setting ACER=1 user decides to generate hard process from AcerMC library. Modeling of ISR/FSR shower, hadroni- 
sation and decays are generated by PYTHIA generator. All steering parameters, relevant for these steps of full event generation 
remain the same as in standard PYTHIA execution. 

By setting ACER=0 user decides to generate standard PYTHIA process. The simple example how to generate tt and ttH 
production process within AcerMC framework is provided in demo_py.f. 

The additional settings ACER=2 and ACER=3, which allow the user to produce event records according to the Les 
Houches standard are given in the Section 7.8. 



'This was however possible in the old PYTHIA 6.1 interface 
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7.6 Details on the interface to HERWIG 6.5 

Interfacing the AcerMC to HERWIG 6.5 is almost identical with the PYTHIA implementation since the interfaced version 
now also complies with the Les Houches standard 1 3 1. The interface routines are written in accordance with the Les Houches 
description; in the HERWIG 6.5 interface the UPINIT routine also has to be called in two steps in order to enable the user to 
change the HERWIG 6.5 default settings and get the correct re-evaluation of the maximal weight (as in the PYTHIA inter- 
face): 

CALL UPINIT(1) 
..user values.. 
CALL UPINIT(2) 

The first UPINIT call is made from the original location in the HWIGUP routine and the second call to UPINIT is 
placed at the end of the HWUINC routine, after all the internal HERWIG settings have been modified. As in the PYTHIA 
implementation all the interface subroutines needed for communication between HERWIG and AcerMC are stored in ac- 
ermc_src/interface/henwig_ac.f. 

One minor additional modification of the original HERWIG 6.5 code was albeit necessary, namely the IMPLICIT NONE 
was commented out in the herwig6500.incfile; this was needed since the AcerMC code is written with the implicit IMPLICIT 
DOUBLE PRECISION(A-H,0-Z). 

In principle the implemented changes should be very easy and transparent for the transfer into new HERWIG releases; an 
example of the use of AcerMC/HERWIG interface is provided in the file demo_hw.f. 

By setting ACER=1 the user decides to generate a hard process from AcerMC library. Modeling of ISR/FSR shower, 
hadronisation and decays are generated by HERWIG generator. All steering parameters, relevant for these steps of full event 
generation remain valid as for the standard HERWIG execution. 

By setting ACER=0 the user decides to generate standard HERWIG processes. The simple example how to generate Zbb 
production process within the AcerMC framework is provided in demo_hw.f. 

The additional settings ACER=2 and ACER=3, which allow the user to produce event records according to the Les 
Houches standard are given in the Section 7.8. 
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The output logs of the run are produced in the directory prod, the acermc.OUt file containing the AcerMC specific in- 
formation and the outputs pythia.out and/or herwig.out listing the outputs of the respective supervising generators. The 
information about the input values of the steering files is stored in run. out in order to facilitate the event generation 'book- 
keeping'. The sample outputs are given in Appendix C. 

Specific information produced by TAUOLA and PHOTOS is stored in respective files tauola.OUt and photOS.OUt. A 
point to stress is that in case the tau decay was restricted the hard process cross-section given in pythia.OUt.herwig.OUt or 
acermc.OUt should be multiplied by a branching ratio as detailed at the end of tauola.OUt, for example: 

< TAUOLA BRANCHING RATIO FOR TAU DECAYS > 



THE TAU DECAYS ARE RESTRICTED TO A: 



LEPTON-HADRON DECAY MODE 



THE PROCESS CROSS-SECTION MUST BE MULTIPLIED BY: 



-> A BRANCHING RATIO = 0.459303E+00 FOR TWO TAUS 



IN THE HARD PROCESS DECAY PRODUCTS!!! 



> TAUOLA BRANCHING RATIO FOR TAU DECAYS < 



7.7 Definition of tlie energy scale 

A few different values of scale used in the evolution of parton density functions as well as the running couplings as{Q^) 
and aQEoiQ^) can be set by the switch ACSET(2) (remember that the factorisation and renormalisation scales are assumed 
to be equal in AcerMC). Note that the correct value of the scale to be used for certain processes is in principle not known; 
what was implemented in AcerMC are the most probable/usual choices on the market; in measurements the best value will 
have to be determined by data analysis. 

• Processes 1,2, 9, 10: 
ACSET(2): (D=1) 

1 - Q= = s 

2 - = E (Pt^ + =<ml> 

3 - = E (Pt')/4 =<pI> 

4-Q2 = (mt + mH/2)^, = 120 GeV/c^ 

• Processes 3 — ^4: 
ACSET(2): (D=1) 

1 - Q' = Mw 

2 - = s*^, where q = b,t 

3 - = + pT^ 

4 - = 0.5 • is'„ + s'^g) + ip^f, where g = b, t 

• Processes 5 — * 8: 

ACSET(2): (D=1) 

1 - Q2 ^ jj^l 

2 - = where q = b,t 
Z-Q'' = sl 

4 - = pT| -I- s*j)/2, where g = b,t 

• Processes 11 —> 14: 
ACSET(2): (D=1) 

1 -Q^ = {2 -ml) 

2- Q' = E(PT'+mJ)/4 

3 - = E (pi'')/2 
4-Q'' = s 

• Processes 15 — > 16: 
ACSET(2): (D=1) 

1 - = s 

2 - = E (Pt" + '"?)/4 = <m^> 

3- g' =E (Pt')/4 =<pi.> 

4- Q2 = (mt + mjj/2)^, = 120 GeV/c^ 
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• Processes 91: 
ACSET(2): (D=1) 
1 - = s 

= (E (p^') + M|)/2 

3-Q' =E (Pt')/2 
4 - q2 = 

« Processes 92 ^ 93: 
ACSET(2): (D=1) 

1- Q^=(2mt)2 =i: (pi,^+m?)/2 

3-Q' =T.iPT^)/2 
=s 

• Processes 94: 
ACSET(2): (D=1) 
^ -Q^ = s 

2- Q^ = {Y: {pir^) + Ml)/2 

3- Q' =E (pi.')/2 
4 - Q2 = 

• Processes 95: 
ACSET(2): (D=1) 

1 -Q^ = (2-m?) 

2-Q' =i:(p^Vm?)/2 

4 - = s 

there are the settings implemented so far. 

7.8 Installation procedure 

The installation requires availability of the CERNLIB fortran library. 

• Ungzip and untar distribution file. 

• In the main directory type make demo_py or make demo.hw. It will compile demo_py.f or demo_hw.f and produce 
the executables demo_py.exe, demo_hw.exe depending on the selected option. The first-time call will also build and 
install all the required libraries; this will not be repeated when the demo_py.f or demo_hw.f are subsequently changed. 

• To execute the programs type make run_py or make run_hw. 

The scripts will change directory to prod and create respective links to data directories there. The execution will also be 
performed there. All input files should be accessible/routed from directory prod, the output files will also be produced 
in that directory. 

7.9 Storing and reading events using the Les Houches accord 

When setting the ACER=2 switch the AcerMC 1.4 is instructed to dump the events generated by the AcerMC library in 
combination with PYTHIA or HERWIG into a pair of output files using the Les Houches format f?']. The files are Ac- 
erMG.pXXX_rYYYYYYY.inparm and Ace rMC_pXXX.rYYYYYYY. events, with XXX denoting the AcerMC process ID 
while YY. is currently the random seed used in AcerMC, to provide a unique identifier in the absence of other choices (run 
number & similar). 

The *.inparm file is basically an information header, containing the relevant run parameters in accordance with the Les 
Houches standard | 3 1, e.g. the cross-section of the process and the number of dumped/stored events; it is written in a 'human- 
readable' format in order to provide the user with the necessary information. The *. events file on the other hand contains the 
actual record of the hard process events produced by the AcerMC. The event dump is obtained by compiling and running the 
provided programs demo_py.exe or demo_hw.exe as detailed in Section lTISl 

Using the ACER=3 switch the thus produced event record can subsequently be read back either into the AcerMC 1.4 
generator to be processed further by PYTHIA or HERWIG, i.e. the ISR/FSR and hadronisation can be added. 

The same event records can also be read into the stand-alone PYTHIA or HERWIG generators, by re-directing the origi- 
nal(dummy) UPINIT routine to the INITACERMG routine and the original UPEVNT routine to USEACERMC routine (In 
the AcerMC 1.4 setup this is done automatically). 
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The files containing the respective routines are named initacermc.f and useacermc.f and are provided in the leshouches 
directory of the AcerMC 1.4 distribution. The reading of a specified set of generated events is achieved by copying the 
*.inparm and *. events files into the running directory and linking/copying the desired *.inparm file to a file named inpar- 
mAcerMC.dat; there is no need to modify the name of the *. events file since its name is aheady stored in the corresponding 
*.inparm file. 

7.10 Interface of tauola to pythia and herwig 

• Interface to PYTHIA : 

The choice of TAUOLA library in PYTHIA sets the t-s to be stable by setting MDCJ(15,1)=0, thus leaving them to 
be treated by TAUOLA procedures. The TAUOLA library is called via the TAUOLA_HEPEVT(IMODE) routine, where 
IM0DE=-1 ,0,1 represent the initialisation, operation and finalisation respectively. In the operation mode the call to 
TAUOLA.HEPEVT(O) is made after the PYEVNT and PYHEPC(1 ) calls which generate the full event and translate and fill 
it into the HEPEVT common block. After TAUOLA has finished a subsequent call PYHEPC(2) is made to translate the new 
event back into the PYTHIA internal structure and a subsequent call to PYEXEC is made so that any undecayed particles 
(e.g. the 7r°) are decayed. It was decided against using the provided PYTHIA interface via the PYTAUD routine since this 
routine is called for each occurrence of a T-lepton separately so no (complex) polarisation options can be included without 
substantial changes to PYTHIA routines. 

• Interface to HERWIG : 

The present version 6.5 of HERWIG contrary to PYTHIA does internal tracking of the polarisations of the r-leptons for 
the more complex built-in processes (e.g. Higgs decays); contrariwise nothing can be done for the external hard processes 
passed to HERWIG for hadronisation. In order to preserve this feature the HWDTAU routine has been modified to perform 
the TAUOLA calls only for internal HERWIG processes which provide the r polarisation information; other r-lepton decays 
are again executed via a call to TAUOLA_HEPEVT. Since HERWIG is using the HEPEVT common block already as the 
internal event record no translation of events is needed. In order to also keep other HERWIG parameters current (e.g. IDHW 
of the particle) a call to a new routine HWHEPC is provided; this new routine also corrects the vertex positions of the r decay 
products and sets the status codes of new particles to the HERWIG recognisable values. Subsequently a call to HWDHAD is 
made in order to decay any undecayed particles. 

7.11 Interface of photos to pythia and herwig 

• Interface to PYTHIA : 

Choosing PHOTOS to provide the QED final state radiation sets the parameter PARJ(90)= 2 • 10'' order to prevent 
PYTHIA to radiate photons off leptons, thus inducing double counting. The parameter is representing the threshold in GeV 
below which leptons do not radiate. Pythia does not contain the routines to handle radiation of hadrons however some caution 
is necessary due to some exceptions/specific decays (e.g. 7) which are generally recognised by PHOTOS itself. 

The initiaUsation, execution and finalisation are done by calls to the new PHOTOS.HEPEVT(IMODE), IMODE=-1,0,1 
subroutine, constructed for this purpose by the AcerMC authors. Since PHOTOS operates on the HEPEVT record the calls 
to PYHEPC routine are again necessary. 

• Interface to HERWIG : 

The present version 6.5 of HERWIG does itself not provide the final state QED radiation of any kind, thus the inclusion of 
PHOTOS is simple and possibly also rather necessary. The call sequence is the same as in PYTHIA , with no need for event 
record conversion. 

7.12 Details of the tauola implementation 

The TAUOLA library is built from the latest distribution source with the Cleo setup option. The native random generator was 
replaced with a link to the random generator of the linked event generator (PYTHIA or HERWIG ) in the same manner as 
done for AcerMC , which decreases the number of random seeds which need to be initialised. 

The TAUOLA_HEPEVT routine is a modification of the original TAUOLA routine provided in the file taufacejetset.f 
in the TAUOLA distribution. The modifications were restricted to allow for the 'overloaded' use of the HEPEVT record by 
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PYTHIA and HERWIG , where respective mother and daughter pointers (which should match) sometimes point to different 
particles (e.g. hard process copies and similar). The original TAUOLA routine already worked when interfaced to PYTHIA 
, albeit requiring the non-default setting of MSTP(1 28)=1 , whereby all the (sometimes useful) links to the hard record were 
lost. The modified version also works with the default PYTHIA MSTP(1 28)=0 setting. 

Additional modifications were made in order to use the parameters set in tauola.card and the call to INIMAS routine was 
replaced by a call to TAU_INIMAS, which is a copy of the former but sets the particle masses to the values of PYTHIA or 
HERWIG defauhs. 

A routine TAUBRS that defines the special decay modes (e.g. LEPTON-HADRON, where one tau decays hadronically 
and the other one leptonically) was written; it operates by modifying the internal GAMPRT array on an event by event basis. 
A related routine TAUBR_PRINT prints the value of the branching ratio into the tauola.OUt file.No changes to the native 
tauola code was necessary (and therefore not made). 

7.13 Details of the photos implementation 

The PHOTOS is also built from the latest distribution source. As in TAUOLA The native random generator was replaced 
with a link to the random generator of the Unked event generator (PYTHIA or HERWIG ) in the same manner as done for 
AcerMC . 

The PHOTOS_HEPEVT routine and the subsequently called PHOTOS.HEPEVT.MAKE are newly written routines 
that inspects the HEPEVT record for charged particles and finds their highest 'mother' particle, which is consequently passed 
to photos routines as the starting point. Photos itself then walks down the branches and performs radiation where possible. A 
bookkeeping of the starting points is made in PHOTOS_HEPEVT in order to prevent multiple invocations of radiation on the 
same particle. 

Original photos code had to be modified due to the 'overloaded' HEPEVT record, since its requirements for matching 
mother-daughter pointers were too strict for either PYTHIA (with external processes and/or l\/ISTP(1 28)=0 setting) or HER- 
WIG . The modification was limited to PHOTOS.MAKE and PHOBOS routines. In addition the tracking of IDHW array 
was added to PHOTOS_MAKE to accommodate the HERWIG event record. A further modification was however necessary 
in the PHOIN routine since in HERWIG the entry JM0HEP(2,I) is not empty but filled with colour flow information, which 
in turn inhibited PHOTOS radiation off participating particles". No modifications of the core (physics) PHOTOS code was 
made. 



"PHOTOS expects the non-zero second 'mother' JM0HEP(2,I) entry only for gg{qq) tt process, which is treated by a set of dedicated routines; 
this in turn clashes with HERWIG 'overloaded' JM0HEP(2,I) entry. 
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8 Outlook and conclusions 

In this paper we presented the AcerMC Monte Carlo Event Generator, based on the hbrary of the matrix-element-based 
generators and interfaces to the universal event generators PYTHIA 6.2 and HERWIG 6.5. The interfaces are based on the 
standard proposed in 1 3 1. 

The presented library fulfills the following goals: 

• It gives a possibility to generate the few Standard Model background processes which were recognised as very dangerous 
for the searches for the New Physics at LHC, and generation of which was either unavailable or not straightforward so 
far. 

• Although the hard process event is generated with matrix-element-based generator, the provided interface allows to 
complete event generation with initial and final state radiation, multiple interaction, hadronisation, fragmentation and 
decays, using implementation of either PYTHIA 6.2 or HERWIG 6.5. 

• These interfaces can be also used for studying systematic differences between PYTHIA 6.2 or HERWIG 6.5 predictions 
for the underlying QCD processes. 

The complete list of the native AcerMC processes implemented so far is: gg, qq — > tibb; qq — ^ W{'^ lv)bb; gg, qq 
Z/^*{^ M)bb; qq lv)tt; gg,qq ^ Z/'y*{^ ££, uiy,bb)tt; gg,qq {Z/W/'y* ^)ttbb; gg,qq tttt; 

gg,qq {tt ~->),f fbf fb ; gg,qq {WWbb —^)ff.ffbb. We plan to extend this not too exhaustive, but very much 
demanded list of processes, in the near future. 

Several improvements of the existing Monte Carlo algorithms/programs have been developed in the process of this work. 
Let us make short list of the most interesting ones: (1) The use of the adapted Kajantie-Byckling enables one to automa- 
tise and modularise the phase space generation of n-particle final states. (2) The additions and extensions to the available 
(multi-channel) phase space algorithms (e.g. Breit-Wigner function with s-dependent width, mass-threshold effects) lead to 
substantial improvement of the unweighting efficiency; Figs.lSi nTil and lTTl illustrate the improvements achieved in the gener- 
ation efficiency. (3) The power of the multi-channel optimisation was enhanced by using the modified ac-VEGAS package. 
We believe that the modification in the VEGAS code represents a very powerful extension of this package; (4)the colour flow 
information has been obtained after some modification of MADGRAPH package. 

Having all these different production processes implemented in the consistent framework, which can be also directly used 
for generating standard processes available in either PYTHIA 6.2 or HERWIG 6.5 Monte Carlo, represents very convenient 
environment for several phenomenological studies dedicated to the LHC physics. Such frame was not available to our knowl- 
edge so far. We hope that it can serve as an interesting example or even a framework. This way some tools for discussing the 
ambiguities due to QCD effects are collected, however the necessary discussion for the appropriate uncertainties is still not 
exhausted. Nevertheless some discussions using this tool can be already found in I, lOJ . LllJ . L12J . L16J . 
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A Feynman Diagrams 

The 38+7 Feynman diagrams contributing to the gg,qq — + ttbb production. Only four flavours are included for incoming 
quarks. Contribution of the incoming b-quarks could be excluded from the calculations thanks to very high suppression 
induced by either the parton density functions and/or CKM matrix elements. 



Diagrams by MadGraph Diagrams by MadGraph 




Diagrams by MadGraph Diagrams by MadGraph 




graph 37 graph 38 graph 7 



Fig. 19: The Feynman diagrams for the processes gg, qq ttbb. 
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The Feynman diagrams contributing to the qq —> W(—> li>)bb and qq —> lv)tt matrix element are just two t- 

chaimel diagrams with fermion exchange and double conversion into an off-shell W boson and a virtual gluon; the W boson 
subsequently decays leptonically into Iv and the gluon splits into a hb pair or ti pair respectively, (c.f. Figure l20t . 



Diagrams by MadGraph 




graph 1 graph 2 



Fig. 20: The Feynman diagrams for the process qq —> Wbb — > e^Ve bb. The same set is used 
for qq — > Wbb e^v,, ti process, with b-quarks replaced by top-quarks. 
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The Feynman diagrams contributing to the gg, qq Z/'y* {ff, vv)hh production are shown in Figure The dominant 
contribution comes from the (2) and (6) configurations for the processes with gg initial state and the double conversion config- 
uration (2), (4) for the ones with qq initial state. The same set of Feynman diagrams is used for the gg, qq — > Z/ 7* (//, vi^)tt 
process, with b-quarks being replaced by the top-quarks. If the Z/'y* bb) decay mode is simulated, it represents only 
subset of the EW production of tibb final state. 



Diagrams by MadGraph Diagrams by MadGraph 




graph 13 graph 14 graph 15 




graph 16 



Fig. 21: The Feynman diagrams for the processes gg, qq — s- Z/"f*bb — > eebb. 
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The complete set of th the Feynman diagrams contributing to the full electro-weak gg, qg —> {Z/W/y* —*)bbtt produc- 
tion mediated by exchange of the Z/W/y* bosons is shown in Fig. 1221 



Diagrams by MadGraph 




Diagrafns by MadGraph 
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The set of the Feynman diagrams contributing to the gg,qq —> if production implemented with different approaches: 
resonant only 2^6 process gg, qq —» {tt —*).ffbffh, complete 2 — > 6 process gg, qq — > {WbWb —>)ffbffb and 2 — > 4 
process gg WbWb — > are shown in Figs. l24l - |27l . 



Diagrams by MadGrapli 9 9 -> ta- vt~ ta+ vt b b- 




Fig. 24: The Feynman diagrams for the processes gg^ qq — > {tt ~^)ffbffb. 



Diagrams by MadGraph 9 9 -> ta+ vt ta- vt- b b- 

, 5 1 



Diagrams by MadGraph 
3 1 . 



> ta+ vtta- vt~ b fa- 




Fig. 25: The Feynman diagrams for the processes gg — ^ {WbWb -—^)ffbffb. 
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Diagrams by MadGraph u u- -> ta+ vt ta- vt- b b~ 
1 ^71 ^'^ 




graph 1 graph 2 graph 3 graph 4 



1 ^71 ^51 ^51 3 




graph 5 graph 6 graph 7 graph 8 



1 ^31, y5 1 ,51 ,7 




graph 9 graph 10 graph 11 graph 12 




graph 13 graph 14 



Fig. 26: r/ie Feynman diagrams for the processes qq {WbWb ~^)ffbffb. 



Diagrams by MadGraph g g -> w- w+ b b- 




Fig. 27: The Feynman diagrams for the processes gg WbWb. 



The 76+14 Feynman diagrams contributing to the gg, qq — * tttt production are shown in Figures |28land|29| Only four 
flavours are included for incoming quarks. Contribution of the incoming b-quarks could be excluded from the calculations 
thanks to very high suppression induced by either the parton density functions and/or CKM matrix elements. 
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Fig. 28: The Feynman diagrams for the processes gg — > tttt. 
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Fig. 29: The Feynman diagrams for the processes q — > tttt. 
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B Example input files 
B.l File run.card 

c 

C STEERING FILE FOR ACERMC (2.0) - ACERMC SETTINGS 

C 

C==== TURN ON FFKEY STEERING FILE (DEBUG) 
LIST 

C THE AcerMC EVENT SETTINGS 

C==== SCALE/ACSET (2) 

C Choose the Q"2 scale for the active AcerMC process. 

C The implemented values differ for various processes, please look into the manual 
C for details 

ACSET2 1 

C==== FERMION/ACSET (3) 

C The flavour of the final state leptons produced in W or Z decays of AcerMC 
C processes . The Pythia/PDG naming convention is used: 

C FERMI0N=1 - only hadronic decays 
C FERMI0N=4 - electron and muon 
C FERMI0N=5 - b-quark decay 
C FERMI0N=11 - electron only 
C FERMI0N=13 - muon only 
C FERMI0N=15 - tau only 

C FERMI0N=12 - neutrinos, all three flavours are generated and 
C the cross-section is calculated accordingly. 
C The setting FERMI0N=5 works only for processes 7 and 8 
ACSET3 13 

C==== Z/GAMMA/ACSET (4) 

C Use the full Z/gamma* propagator in AcerMC processes 5-8. 

C ZGAMMA=0 - only Z propagator 

C ZGAMMA=1 - full Z/gamma* propagator 

ACSET4 1 

C==== Z/G CUT/ACSET(5) 

C Cutoff value on the invariant mass m_Z/gamma* in GeV when ZGAMMA=1 . 

C Note that the provided data files exist only for values of 

C ZGCUT=10, 30, 60, 120, 300 GeV which should satisfy most 

C users. In case a different value is set the user has also to provide 

C the user data files for the run. 

ACSET5 60.0 

C THE AcerMC ADVANCED SWITCHES 

C==== ALPHA_S/ACSET (8) 

C Use the alpha_s provided by the linked generator (Pythia/Herwig) or the 
C one provided by AcerMC 
herwig . pdf 

C ALPHAS-0 - use the linked generator's alpha_s 
C ALPHAS=1 - use the AcerMC s alpha_s (one loop calculation) 
C ALPHAS-2 - use the AcerMC s alpha_s (three loop calculation) 
ACSET8 

C==== ALPHA_S (M_Z~2) /ACSET (9) 

C Specify the value of alpha_QCD (M_Z"2) for AcerMC alpha_s calculation 
C ALPHASMZ=-1 - the value is taken from PDFLIB 
C ALPHASMZ>0 - the provided value is taken 
ACSET9 -1. 
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C==== ALPHA_EM/ACSET (10) 

C Use the alpha_QED provided by the linked generator (Pythia/Herwig) or the 

C one provided by AcerMC 

C ALPHAEM-0 - use the linked generator' s alpha_QED 
C ALPHAEM=1 - use the AcerMC s alpha_QED 
ACSETIO 

C==== ALPHA_EM(0) /ACSET (11) 

C Specify the value of alpha_QED(0) for AcerMC alpha_QED calculation 
C ALPHAEM0=-1 - the default AcerMC value is used 
C ALPHAEMO>0 - the provided value is taken 

ACSETll -1. 

C==== TOP S-L/ACSET (12) 

C Specify the decay mode of WW pair produced by external top decays in AcerMC 
C processes 1,2,4,7,8,9,92 and 93: 

C TOPDEC=0 - both W bosons decay according to Pythia/Herwig switches 

C T0PDEC=1 - one W decays into electron + nu and the other one hadronically 

C T0PDEC=2 - one W decays into muon + nu and the other one hadronically 

C T0PDEC^3 - one W decays into tau + nu and the other one hadronically 

C T0PDEC=4 - one W decays into el or mu + nu and the other one hadronically 

C TOPDEC^S - one W decays into el or mu + nu and the other one hadronically, the 

C W decaying leptonically has the same charge as the primary W; the decay 

C mode makes sense only for AcerMC processes 4! 

C When TOPDEOO the output cross-section is ALREADY MULTIPLIED by the corresponding 
C branching ratio (s) ! (Courtesy of AcerMC authors) 

ACSET12 

C==== BOSON PAIR DECAYS/ACSET (13) 

C Specify the decay mode of boson pairs inside AcerMC processes: 
C 

C BOSDEC^O - both bosons decay in all possible modes 

C B0SDEC=1 - one boson decays into electron + nu and the other one hadronically 

C B0SDEC^2 - one boson decays into muon + nu and the other one hadronically 

C B0SDEC^3 - one boson decays into tau + nu and the other one hadronically 

C B0SDEC=4 - one boson decays into el or mu + nu and the other one hadronically 

C B0SDEC=11 - both bosons decay into el + nu 

C B0SDEC=13 - both bosons decay into muon + nu 

C B0SDEC=15 - both bosons decay into tau + nu 

C BOSDEC=20 - both bosons decay hadronically 

C The output cross-section is ALREADY MULTIPLIED by the corresponding 
C branching ratio (s) ! 

ACSET13 13 

C THE AcerMC TRAINING SETUP AND UNWEIGHTING TREATMENT 

C==== MODE/ACSET (50) 

C Specify the AcerMC training mode: 

C MODE-0 - normal run, generate unweighted events 
C MODE^l - perform multi-channel optimisation . 
C MODE-2 - perform VEGAS grid training. 

C M0DE=3 - perform VEGAS grid training as M0DE=2 but does this by upd 

C M0DE=-1 - all flat 

C M0DE=-2 - only VEGAS flat 

ating a provided grid 

ACSET50 

C==== USER/ACSET (52) 

C Use the data files provided by user 
C USER=0 - no, use internal files 

C USER=1 - use the user' s multi-channel optimisation and VEGAS grid files 

C USER=2 - use the default multi-channel optimisation and user' s VEGAS grid files 

C USER=3 - use the default multi-channel optimisation and VEGAS grid files, read the user 

C maximal weight file . 
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ACSET52 

C==== MAXFIND/ACSET (53) 

C Search for the maximum weight needed for event unweighting 
C MAXFIND-0 - no, use the provided file for max. weights 

C MAXFIND=1 - use the provided file for max. weights, re-calculate the max. weights using 
C the stored 100 highest events 

C MAXFIND=2 - perform the search and give the wtmax file, equivalent to generation of 
C weighted events 

ACSET53 

C==== EPSILON/ACSET (54) 

C Use the epsilon maximal weight or the overall maximal weight found in training (see the 

C manual for the difference) 

C EPSILON=0 - use the epsilon max. weight 

C EPSIL0N=1 - use the overall maximal weight 

ACSET54 

C==== NQCD/ACSET (56) 

C Use the naive qcd correction for width calculations 

C (see the manual for details) 

C NQCD=1 - use the naive QCD corrections 

C NQCD=0 - don't use naive QCD corrections 

ACSET5 5 1 

Q 

END 

B.2 File acermccard 

c 

C STEERING FILE FOR ACERMC (2.0) - ACERMC SETTINGS 

C 

C==== TURN ON FFKEY STEERING FILE (DEBUG) 
LIST 

C THE AcerMC EVENT SETTINGS 

C==== SCALE/ACSET (2) 

C Choose the Q"2 scale for the active AcerMC process. 

C The implemented values differ for various processes, please look into the manual 
C for details 

ACSET2 1 

C==== FERMION/ACSET (3) 

C The flavour of the final state leptons produced in w or z decays of AcerMC 

C processes . The Pythia/PDG naming convention is used: 

C FERMI0N=1 - only hadronic decays 

C FERMI0N=4 - electron and muon 

C FERMI0N=5 - b-quark decay 

C FERMI0N=11 - electron only 

C FERMI0N=13 - muon only 

C FERMI0N=15 - tau only 

C FERMI0N=12 - neutrinos, all three flavours are generated and 
C the cross-section is calculated accordingly. 
C The setting FERMI0N=5 works only for processes 7 and 8 
ACSET3 13 

C==== Z/GAMMA/ACSET (4) 

C Use the full z/gamma* propagator in AcerMC processes 5-8. 

C ZGAMMA=0 - only Z propagator 

C ZGAMMA=1 - full Z/gamma* propagator 
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ACSET4 1 

C-— Z/G CUT/ACSET{5) 

C Cutoff value on the invariant mass m_Z/ gamma* in GeV when ZGAMMA=1 . 
C Note that the provided data files exist only for values of 
C ZGCUT=10, 30, 60, 120, 300 GeV which should satisfy most 

C users. In case a different value is set the user has also to provide 
C the user data files for the run. 

ACSET5 60.0 

C THE AcerMC ADVANCED SWITCHES 

C==== ALPHA_S/ACSET (8) 

C Use the alpha_s provided by the linked generator (Pythia/Herwig) or the 
C one provided by AcerMC 
herwig . pdf 

C ALPHAS=0 - use the linked generator' s alpha_s 
C ALP HAS ^1 - use the AcerMC s alpha_s (one loop calculation} 
C ALPHAS=2 - use the AcerMC s alpha_s (three loop calculation) 
ACSET8 

C==== ALPHA_S (M_Z"2) /ACSET (9) 

C Specify the value of alpha_QCD (M_Z " 2 ) for AcerMC alpha_s calculation 
C ALPHASMZ— 1 the value is taken from PDFLIB 
C ALPHASMZ>0 - the provided value is taken 
ACSET9 -1. 

C==== ALPHA_EM/ACSET (10) 

C Use the alpha_QED provided by the linked generator (Pythia/Herwig) or the 
C one provided by AcerMC 

C ALPHAEM-0 - use the linked generator's alpha_QED 
C ALPHAEM-1 - use the AcerMC s alpha_QED 
ACSETIO 

C==== ALPHA_EM(0) /ACSET (11) 

C Specify the value of alpha_QED (0) for AcerMC alpha_QED calculation 
C ALPHAEM0=-1 - the default AcerMC value is used 
C ALPHAEMO>0 - the provided value is taken 

ACSETll -1. 

C==== TOP S-L/ACSET {12} 

C Specify the decay mode of WW pair produced by external top decays in AcerMC 
C processes 1,2,4,7,8,9,92 and 93: 

C TOPDEC=0 - both W bosons decay according to Pythia/Herwig switches 

C T0PDEC=1 - one W decays into electron + nu and the other one hadronically 

C T0PDEC^2 - one W decays into muon + nu and the other one hadronically 

C TOPDEC^S - one W decays into tau + nu and the other one hadronically 

C T0PDEC=4 - one W decays into el or mu + nu and the other one hadronically 

C T0PDEC=5 - one W decays into el or mu + nu and the other one hadronically, the 

C W decaying leptonically has the same charge as the primary W; the decay 

C mode makes sense only for AcerMC processes 4! 

C When TOPDEOO the output cross-section is ALREADY MULTIPLIED by the corresponding 
C branching ratio (s) ! (Courtesy of AcerMC authors) 

ACSET12 

C — - BOSON PAIR DECAYS/ACSET (13) 

C Specify the decay mode of boson pairs inside AcerMC processes: 
C 

C BOSDEC=0 - both bosons decay in all possible modes 

C B0SDEC=1 - one boson decays into electron + nu and the other one hadronically 
C B0SDEC=2 - one boson decays into muon + nu and the other one hadronically 
C B0SDEC=3 - one boson decays into tau + nu and the other one hadronically 
C B0SDEC=4 - one boson decays into el or mu + nu and the other one hadronically 
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C B0SDEC=11 - both bosons decay into el + nu 
C B0SDEC=13 - both bosons decay into muon + nu 

C BOSDEC^IS - both bosons decay into tau + nu 
C BOSDEC=20 - both bosons decay hadronically 

C The output cross-section is ALREADY MULTIPLIED by the corresponding 
C branching ratio (s) ! 

ACSET13 13 

C THE AcerMC TRAINING SETUP AND UNWEIGHTING TREATMENT 

MODE/ACSET (50} 
C Specify the AcerMC training mode: 

C MODE=0 - normal run, generate unweighted events 
C MODE^l - perform multi-channel optimisation . 
C M0DE=2 - perform VEGAS grid training. 

C M0DE=3 - perform VEGAS grid training as M0DE=2 but does this by upd 

C M0DE=-1 - all flat 

C M0DE=-2 - only VEGAS flat 

ating a provided grid 

ACSET50 

C==== USER/ACSET (52) 

C Use the data files provided by user 
C USER=0 - no, use internal files 

C USER=1 - use the user' s multi-channel optimisation and VEGAS grid files 

C USER=2 - use the default multi-channel optimisation and user' s VEGAS grid files 

C USER=3 - use the default multi-channel optimisation and VEGAS grid files, read the user 

C maximal weight file . 

ACSET52 

MAXFIND/ACSET (53) 
C Search for the maximum weight needed for event unweighting 
C MAXFIND=0 - no, use the provided file for max . weights 

C MAXFIND=1 - use the provided file for max. weights, re-calculate the max. weights using 
C the stored 100 highest events 

C MAXFIND=2 - perform the search and give the wtmax file, equivalent to generation of 
C weighted events 

ACSET53 

EPSILON/ACSET (54) 

C Use the epsilon maximal weight or the overall maximal weight found in training (see the 

C manual for the difference) 

C EPSILON=0 - use the epsilon max. weight 

C EPSIL0N=1 - use the overall maximal weight 

ACSET54 

C==== NQCD/ACSET (56) 

C Use the naive qcd correction for width calculations 

C (see the manual for details) 

C NQCD=1 - use the naive QCD corrections 

C NQCD=0 - don't use naive QCD corrections 

ACSET5 6 1 

C 

END 



B.3 File tauola.card 



C STEERING FILE FOR TAUOLA & ACERMC (1.4) 
C 
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C==== TURN ON FFKEY STEERING FILE (DEBUG) 
LIST 

C==== POLAR 

C Polarisation switch for tau decays 
C POLAR=0 - switch polarisation off 
C POLAR-1 - switch polarisation on 
POLAR 1 

C==== RADCOR 

C Order (alpha) radiative corrections for tau decays 

C RADCOR=0 - switch corrections off 
C RADCOR-1 - switch corrections on 

RADCOR 1 
C==== PHOX 

C Radiative cutoff used in tau decays 

C PHOX=0.01 - default value by TAUOLA authors 

PHOX 0.01 
C==== DMODE 

C Tau and tau pair decay mode 

C DMODE=0 - all decay modes allowed 

C DM0DE=1 - (LEPTON-LEPTON) : only leptonic decay modes 
C DM0DE=2 - (HADRON-HADRON) : only hadronic decay modes 

C DM0DE=3 - (LEPTON-HADRON) : one tau decays into leptons and the other one into hadrons 
C DM0DE=4 - (TAU->PI NU) : taus are restricted to decay to a pion and neutrino 

DMODE 2 

C==== JAK1/JAK2 

C Decay modes of taus according to charge, the list is taken from TAUOLA output 
C The listing gives only Tau- modes, the Tau+ are charge conjugate, neutrinos 
C are omitted. 



c 


JAKl/2 




1 




TAU- 


--> 


E- 






c 


JAKl/2 




2 




TAU- 


— > 


MU- 






c 


JAKl/2 




3 




TAU- 


— > 


PI- 






c 


JAKl/2 




4 




TAU- 


— > 


PI-i 


, PIO 




c 


JAKl/2 




5 




TAU- 


— > 


Al- 


(two 


subch) 


c 


JAKl/2 




6 




TAU- 


--> 


K- 






c 


JAKl/2 




7 




TAU- 


--> 


K*- 


(two 


subch) 


c 


JAKl/2 




8 




TAU- 


--> 


2PI- 


PI 


0, PI + 


c 


JAKl/2 




9 




TAU- 


--> 


3PI0, 


PI- 


c 


JAKl/2 




10 




TAU- 


--> 


2PI- 


PI+, 2PI0 


c 


JAKl/2 




11 




TAU- 


— > 


3PI- 


-, 2PI + 


c 


JAKl/2 




12 




TAU- 


— > 


3PI- 


-, 2PI + , PIO 


c 


JAKl/2 




13 




TAU- 


— > 


2PI- 


-, PI+, 3PI0 


c 


JAKl/2 




14 




TAU- 


— > 


K-, 


PI-, 


K+ 


c 


JAKl/2 




15 




TAU- 


--> 


KO, 


PI-, 


KOB 


c 


JAKl/2 




16 




TAU- 


--> 


K-, 


KO, 


PIO 


c 


JAKl/2 




17 




TAU- 


--> 


PIO 


PIO 


K- 


c 


JAKl/2 




18 




TAU- 


--> 


K- 


PI- 


PI + 


c 


JAKl/2 




19 




TAU- 


--> 


PI- 


KOB 


PIO 


c 


JAKl/2 




20 




TAU- 


— > 


ETA 


PI- 


PIO 


c 


JAKl/2 




21 




TAU- 


— > 


PI- 


PIO 


GAM 


c 


JAKl/2 




22 




TAU- 


— > 


K- 


KO 





C==== DECAY MODE OF TAU+ 
JAKl 13 

C==== DECAY MODE OF TAU- 
JAK2 13 

Q 

END 
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jPr4___FJi^_RlM>tos.card 



C STEERING FILE FOR PHOTOS S ACERMC (1.4) 

C 

C==== TURN ON FFKEY STEERING FILE (DEBUG) 
LIST 

C==== PMODE 

C Radiation mode of photos 

C PMODE-1 - enable radiation of photons for leptons and hadrons 
C PM0DE=2 - enable radiation of photons for taus only 
C PM0DE=3 - enable radiation of photons for leptons only 
PMODE 1 



C==== XPHCUT 

c Infrared cutoff for photon radiation 

C XPHCUT=0.01 - default value by PHOTOS authors 

XPHCUT 0.01 



C==== ALPHA 

C Alpha (QED) value 

C ALPHA < - leave default (0.00729735039) 
ALPHA -1. 

C==== INTERF 

C Photon interference weight switch 
C INTERF = 1 - interference is switched on 
C INTERF = - interference is switched off 
INTERF 1 

C==== I SEC 

C Double bremsstrahlung switch 

C ISEC=1 - double bremsstrahlung is switched on 
C ISEC=0 - double bremsstrahlung is switched off 
ISEC 1 

C==== IFTOP 

C Switch for gg(qq)->tt~ process radiation 
C IFT0P=1 - the procedure is is switched on 
C IFTOP=0 - the procedure is is switched off 



IFTOP 
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C Example output files 
C.l File acermcout 



. j%3] :, 



_|xx%xxxx%%%%%+ " 
: ~] %xxxx] xx%x_, +_x_% ^ 
_| |x| |xxl+]]]+]]]]x|xxx] ' 

-+%%xxxx] ]]]+]]++]+] x] I >- 

-+x]]]] 1+]+]+]++]++] l+l+l ]+ 
- - !=++]+]=++++++=!++]=+- 




40000L, lOOOOi jOOOS .aOOOOOL#0 

— ?##0L .aaaa aa .aaaa; aaaa, _aaa, -OOOA _0001- _dO ! ' -400 

dO 40, _W0#V9N0#S dO#V*N#0, 0##0LW0@4#@' 00 j#; JO [ 01 dO' 40 

JOl -#W #0' ?#W ##~ -#0; j##9 00 4#|01|01 00 

_00yyyW0L :0f :000###00001 j#l 00 ?#0@~|01 #0 

##!!!! !#0; -OA _ -OA j#l 00 HH< |01 jOL 



ajOAa 4#Aaa_aj#0' ?OLaa_aaaOL aaJOLaaa, 



_aOOaa _ajOLa *OAaa_aad 



HHHRHl HHHHH 



'9##009! 



'9NW00@ ! 



HHHHRHHHl :HHHHH 7HHHRH ?!##00P! 



AcerMC 2.0 (April 2004), B. P. Kersevan, E. Richter-Was 



< ACTIVATED PROCESSES >■ 



1) 


g 


+ 


g 


-> 


t t" b b" (MG) 


ON 


2) 


q 


+ 


q" 


-> 


t t" b b" (MG) 


OFF 


3) 


q 


+ 


q" 


-> 


(W->) 1 nu_l b b" (MG) 


OFF 


4) 


q 


+ 


q" 


-> 


(W->) 1 nu_l t t" (MG) 


OFF 


5) 


g 


+ 


g 


-> 


(Z0->) 1 1^ b b~ (MG) 


OFF 


6) 


q 


+ 


q~ 


-> 


(Z0->) 1 1~ b b" (MG) 


OFF 


7) 


g 


+ 


g 


-> 


(Z0->) f f t t~ (MG) 


OFF 


8) 


q 


+ 


q~ 


-> 


(Z0->) f t t~ (MG) 


OFF 


9) 


g 


+ 


g 


-> 


{ZO/W/gamma->) t t" b b" 


(MG) OFF 


10) 


q + q- 


-> 


(ZO/W/gamina->) t t" b b" 


(MG) OFF 


11) 


g 


+ 


g 


-> 


t t" off-shell (MG) 


OFF 


12) 


q 


+ 


q" 


-> 


t t" off-shell (MG) 


OFF 


13) 


g 


+ 


g 


-> 


(W W b b" ->) (MG) 


OFF 


14) 


q 


+ 


q" 


-> 


(W W b b~ ->) (MG) 


OFF 


15) 


g 


+ 


g 


-> 


t t~ t t~ (MG) 


OFF 


16) 


q 


+ 


q~ 


-> 


t t~ t t" (MG) 


OFF 


91) 


q 


+ 


q~ 


-> 


(ZO/gamma->) 1 1" (MG) 


OFF 


92) 


g 


+ 


g 


-> 


t t" (MG) 


OFF 


93) 


q 


+ 


q" 


-> 


t t" (MG) 


OFF 


94) 


q 


+ 


q" 


-> 


(W->) 1 nu_l (MG) 


OFF 


95) 


g 


+ 


g 


-> 


W+ W- b b" (MG) 


OFF 
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-< ACERMC SETTINGS >- 



C.M.S ENERGY = 
SCALE CHOICE = 
ACERMC ALPHA_QCD 
ALPHA_QCD (M_Z) 
ACERMC ALPHA_QED 
ALPHA_QED (0) 
TOP->W S-L DECAY 
BOSON PAIR DECAY 



14000.00 
1 


-1.000000 


-1.000000 

13 



[ACSET (1) ] 
[ACSET (2) ] 
[ACSET (8) ] 
[ACSET (9) ] 
[ACSET (10) ] 
[ACSET (11) ] 
[ACSET (12) ] 
[ACSET (13) ] 



OPTIMIZATION 
OPTIM. STEPS 
USER FILES 
MAX . SEARCH 
EPSILON CUTOFF 
NAIVE QCD 





100 



1 



[ACSET (50) ] 
[ACSET (51) ] 
[ACSET (52) ] 
[ACSET (53) ] 
[ACSET (54) ] 
[ACSET (5 5) ] 



READ MAXIMUM WEIGHT (MB) 
READ EPSILON WEIGHT (MB) 



0.236905E-07 
0.108907E-07 



SET MAXIMUM WEIGHT (MB) 
SET WEIGHT CORRECTION 



0.108907E-07 
O.lOOOOOE+01 



-< FINALIZATION FOR PROCESS: 1 >- 
< WEIGHT SURVEY > 



-< TOTAL STATISTICS >- 



CROSS-SECTION ESTIMATE 



VARIANCE ESTIMATE 



4-/- 



0.411408E+01 PB 
0.153340E+00 PB 
0.235132E-01 PB'2 
0.159711E-02 PB-2 



MAXIMUM WEIGHT 
NO. WEIGHTS NE 

NO. WEIGHTS EQ = 
NO. WEIGHTS LT = 
EFFICIENCY FOR ALL WEIGHTS 
EFFICIENCY FOR NONZERO WEIGHTS 



0.975277E-08 

252 






42.141 % 
42.141 % 



NO. WEIGHTS ABOVE EPSILON-CUT 



■> WEIGHT SURVEY <■ 
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* * 3^"!^:^ ■fiyjf'^Ht^jQJltt ******************************************************* 



Welcome to the Lund Monte Carlo! 

PPP Y Y TTTTT H H III A 

PPYY T H HI AA 

PPP Y T HHHHH I AAAAA 

P Y THHIAA 

P Y T H H III A A 

This is PYTHIA version 6.214 
Last date of change: 22 Jan 2003 

Now is 6 Feb 2003 at 18:39:46 

Disclaimer: this program comes 
without any guarantees . Beware 
of errors and use common sense 
when interpreting results . 

Copyright T. Sjostrand (2003) 



** An archive of program versions and documentation is found on the web: 

** http : //www . thep . lu.se/~torb jo rn/Pythia. html 

-k -k 

** When you cite this program, currently the official reference is 

** T. Sjostrand, P. Eden, C. Friberg, L. Lonnblad, G. Miu, S. Mrenna and 

** E. Norrbin, Computer Physics Commun. 135 (2001) 238. 

** The large manual is 

** T. Sjostrand, L. Lonnblad and S. Mrenna, LU TP 01-21 [hep-ph/01082 64 ] . 

** Also remember that the program, to a large extent, represents original 

** physics research . Other publications of special relevance to your 

** studies may therefore deserve separate mention. 

** Main author : Torb jorn S jostrand; Department of Theoretical Physics 2 , 

** Lund University, Solvegatan 14A, S-223 52 Lund, Sweden; 

** phone: + 46 - 46 - 222 48 16; e-mail: torbjorn@thep.lu.se 

** Author : Leif Lonnblad; Department of Theoretical Physics 2 , 

** Lund University, Solvegatan 14A, S-223 62 Lund, Sweden; 

** phone: + 46 - 46 - 222 77 80; e-mail: leif@thep.lu.se 

** Author: Stephen Mrenna; Computing Division, Simulations Group, 

** Fermi National Accelerator Laboratory, MS 234, Batavia, IL 60510, USA; 

** phone : + 1 - 630 - 840 - 2556; e-mail : mrenna@fnal . gov 

** Author : Peter Skands; Department of Theoretical Physics 2 , 

** Lund University, Solvegatan 14A, S-223 62 Lund, Sweden; 

** phone: + 46 - 46 - 222 31 92; e-mail: zeiler@thep.lu.se 



PYINIT: initialization of PYTHIA routines 



I I 

I PYTHIA will be initialized for p+ on p+ user configuration I 

I with 7000.000 GeV on 7000.000 GeV beam energies I 

I I 

I corresponding to 14000.000 GeV center-of-mass energy I 

I I 



******** PYMAXI : summary of differential cross-section maximum search ******** 



I I SUB Subprocess name 
I 



I I 
I Maximum value I 
I I 



User process 601 



1 . 9366E-08 



B. P. Kersevan,E. Richter-Was: 



AcerMC 2.0 Event Generator 



79 



********************** PYINIT: initialization completed ********************** 
Event listing of user process at input (simplified) 



I 


1ST 


ID 


Mothe 




Colours 


P- 


_x 


P- 


.y 


P- 


.z 


E 


m 


1 


-1 


21 








501 


503 





000 





000 


332 


784 


332.784 


0.000 


2 


-1 


21 








502 


504 





000 





000 


-614 


582 


614 .582 


0.000 


3 


1 


6 


1 


2 


501 





81 


702 


20 


957 


-203 


328 


281.214 


175.000 


4 


1 


-6 


1 


2 





504 


28 


925 


29 


174 


-310 


224 


358.541 


175.000 


5 


1 


5 


1 


2 


502 





-52 


937 


-77 


151 


33 


793 


99.598 


4.800 


6 


1 


-5 


1 


2 





503 


-57 


690 


27 


020 


197 


950 


208.013 


4.800 



Event listing (summary) 



I 


particle/ jet 


KS 


KF 


orig 


P- 


_x 


P- 


-Y 


P- 


.z 




E 


m 


1 


!p+! 


21 


2212 








.000 





000 


7000 


000 


7000 


.000 





. 938 


2 


!p+! 


21 


2212 








000 





000- 


7000 


000 


7000 


. 000 





. 938 


3 


!g! 


21 


21 


1 


-0 


290 





974 


332 


784 


332 


.785 





. 000 


4 


!g! 


21 


21 


2 


-0 


177 





335 


-614 


582 


614 


.582 





. 000 


5 


!g! 


21 


21 


3 


-0 


290 





974 


332. 


784 


332 


.785 





. 000 


5 


!g! 


21 


21 


4 


-0 


.177 





335 


-514 


582 


614 


.582 





.000 


7 


!t ! 


21 


6 





81 


.598 


21 


204 


-203 


329 


281 


.203 


175 


.000 


8 


!tbar! 


21 


-6 





28 


808 


29 


427 


-310 


250 


358 


. 575 


175 


. 000 


9 


!b! 


21 


5 





-53 


004 


-76 


938 


33 


870 


99 


.495 


4 


. 800 


10 


!bbar! 


21 


-5 





-57 


868 


27 


617 


197 


911 


208 


.094 


4 


. 800 



11 


t 


A 


2 


6 


7 


81 


598 


21 


204 


-203 


329 


281 


203 


175 


000 


12 


ud_0 


V 


1 


2101 


1 


-0 


058 


-0 


344 


4727 


872 


4727 


872 





579 


13 


tbar 


A 


2 


-6 


8 


28 


808 


29 


427 


-310 


250 


358 


575 


175 


000 


14 


u 


V 


1 


2 


2 





233 


-0 


275 


-19 


002 


19 


009 





330 


15 


b 


A 


2 


5 


9 


-53 


004 


-75 


938 


33 


870 


99 


495 


4 


800 


15 


ud_l 


V 


1 


2103 


2 


-0 


056 


-0 


061- 


6366 


410 


6366 


410 





771 


17 


bbar 


A 


2 


-5 


10 


-57 


868 


27 


617 


197 


911 


208 


094 


4 


800 


18 


u 


V 


1 


2 


1 





348 


-0 


630 


1939 


344 


1939 


344 





330 



sum: 2.00 0.00 0.00 0.01 14000.00 14000.00 



********* PYSTAT: Statistics on Number of Events and Cross-sections ********* 



I I II 

I Subprocess I Number of points I Sigma I 

I I II 

I 1 1 (mb) I 

I I II 

I N:o Type I Generated Tried I I 

I I II 

I I II 

I All included subprocesses I 100 491 I 3.913E-09 I 

I 4 User process 501 I 100 491 I 3.913E-09 I 

I I II 
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C.3 File herwig.out 



HERWIG 6.500 



16 Oct 2002 



Please reference: G. Marchesini, B.R. Webber, 
G.Abbiendi, I.G.Knowles, M.H. Seymour & L. Stance 

Computer Physics Communications 67 (1992) 465 
and 

G.Corcella, I.G.Knowles, G.Marchesini, S.Moretti, 
K.Odagiri, P.Richardson, M.H.Seymour & B.R.Webber, 
JHEP 0101 (2001) 010 

INPUT CONDITIONS FOR THIS RUN 



BEAM 1 (P ) MOM. = 

BEAM 2 (P ) MOM. = 

PROCESS CODE (IPROC) 
NUMBER OF FLAVOURS 
STRUCTURE FUNCTION SET = 
AZIM SPIN CORRELATIONS = 
AZIM SOFT CORRELATIONS = 
QCD LAMBDA (GEV) 
DOWN QUARK MASS 

UP QUARK MASS 

STRANGE QUARK MASS 
CHARMED QUARK MASS 
BOTTOM QUARK MASS 
TOP QUARK MASS 

GLUON EFFECTIVE MASS 
EXTRA SHOWER CUTOFF (Q)= 
EXTRA SHOWER CUTOFF (G)= 
PHOTON SHOWER CUTOFF 
CLUSTER MASS PARAMETER = 
SPACELIKE EVOLN CUTOFF = 
INTRINSIC P-TRAN (RMS) = 
DECAY SPIN CORRELATIONS^ 
SUSY THREE BODY ME 
SUSY FOUR BODY ME 
MIN MTM FRAC FOR ISR =1 , 
1-MAX MTM FRAC FOR ISR =1 , 



7000.00 
7000.00 
-601 



T 
T 

0.1800 
0.3200 
0.3200 

. 5000 
1 . 5500 
4 .8000 



175.0000 
0.7500 
0.4800 
0.1000 
0.4000 
3.3500 
2.5000 
0.0000 
T 
T 
F 

OOOOE-04 
OOOOE-06 



NO EVENTS WILL BE WRITTEN TO DISK 



B_d: Delt-M/Gam =0.7000 Delt-Gam/2*Gam =0.0000 
B_s: Delt-M/Gam = 10.00 Delt-Gam/2*Gam =0.2000 



PDFLIB USED FOR BEAM 1: SET 55 OF CTEQ 
PDFLIB USED FOR BEAM 2: SET 55 OF CTEQ 



Checking consistency of particle properties 



Checking consistency of decay tables 



CHECKING SUSY DECAY MATRIX ELEMENTS 



INPUT EVT WEIGHT 
INPUT MAX WEIGHT 



l.OOOOE+00 
O.OOOOE+00 



EVENT 
SEEDS : 



1: 7000.00 GEV/C P 
725123 S 1450246 



ON 7000.00 GEV/C P 
STATUS: 4 ERROR: 



PROCESS: -601 
WEIGHT: l.OOOOE+00 



-INITIAL STATE- 



IHEP 


ID 


IDPDG 


1ST 


MO I 


MO 2 


DAI 


DA2 


P-X 


P-Y 


P-Z ENERGY MASS 




v-x 


1 


P 


2212 


101 














.00 


0.00 


7000.0 7000.0 0.94 





. OOOE + 00 


2 


P 


2212 


102 














0.00 


0.00- 


-7000.0 7000.0 0.94 


0, 


.OOOE+00 


3 


CMF 





103 


1 


2 








0.00 


0.00 


0.014000.014000.00 


0, 


.OOOE+00 



HARD SUBPROCESS 
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IHEP 


ID 


IDPDG 


1ST 


MOl 


MO 2 


DAI 


DA2 


P-X 


P-Y 


P-Z 


ENERGY 


MASS 




v-x 


4 


GLUON 


21 


121 


6 


10 


11 


7 


.00 


0.00 


204 .5 


204 .5 





00 





OOOE+00 


5 


GLUON 


21 


122 


6 


8 


17 


9 


.00 


0.00 


-266.7 


266 .7 





00 





OOOE+00 


6 


HARD 





120 


4 


5 


7 


10 


-33 . 00 


-38 .11 


-62.2 


474 . 


467 


16 





OOOE+00 


7 


TQRK 


6 


123 


6 


4 


26 


8 


-5.57 


54.25 


82.7 


201.1 


175 


00 





OOOE+00 


8 


TEAR 


-6 


124 


6 


7 


28 


5 


-10.74 


-81.04 


-127.9 


231.7 


175 


00 





OOOE+00 


9 


BQRK 


5 


124 


5 


5 


30 


10 


12.62 


23.88 


-11.1 


29.6 


4 


80 





OOOE+00 


10 


BEAR 


-5 


124 


6 


9 


32 


4 


3.70 


2.91 


-5.8 


8.9 


4 


80 





OOOE+00 



OUTPUT ON LES HOUCHES EVENTS 



PROG CODE XSECT(pb) 
-601 0.29506E+01 



XERR (pb) 
0.10898E+00 



OUTPUT ON ELEMENTARY PROCESS 



N.B. NEGATIVE WEIGHTS NOT ALLOWED 



NUMBER OF EVENTS 
NUMBER OF WEIGHTS 
MEAN VALUE OF WGT 
RMS SPREAD IN WGT 
ACTUAL MAX WEIGHT 
ASSUMED MAX WEIGHT 

PROCESS CODE IPROC 

CROSS SECTION (PB) 
ERROR IN C-S (PB) 
EFFICIENCY PERCENT 



100 

490 

2 . 9506E-03 
O.OOOOE+00 
1.5177E-02 
1.5177E-02 



-601 



2 . 951 
.1090 
19.44 



Max wgt (nb) No. of events 
0.15177E-01 100 



